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Abstract 

We address a model for adhesive unilateral frictionless Signorini-type contact between bodies of heat- 
conductive viscoelastic material, in the linear Kelvin- Voigt rheology, undergoing thermal expansion. The 
flow-rule for debonding the adhesion is considered rate-independent and unidirectional, and a thcrmody- 
namically consistent model is derived and analysed as far as the existence of a weak solution is concerned. 
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1. Introduction 

We are interested in the modelling of elastic bodies glued together by an adhesive, which can undergo 
an inelastic process of so-called delamination (sometimes also called debonding). "Microscopically" speak- 
ing, some macromolecules in the adhesive may break upon loading and we assume that they can never 
be glued back, i.e., no "healing" is possible. This makes the process unidirectional; sometimes it is also 
referred to as irreversible, although this adjective has an alternative thermodynamical meaning as dissi- 
pative in general. On the glued surface, we consider the delamination process as rate- independent and, 
in the bulk, we also consider rate-dependent inertial, viscous-like, and thermal- expansion effects. More- 
over, we confine ourselves to small strains and, just for the sake of notational simplicity, we restrict the 
analysis to the case of two bodies f2+ and Vl- glued together along the contact surface Tq- The material 
in the bulk is considered as heat conductive, and thus the system is completed by the nonlinear heat 
equation in a thermodynamically consistent way. The contact surface is considered infinitesimally thin, 
so that the thermal capacity of the adhesive is naturally neglected. The coupling of the mechanical and 
thermal effects thus results from thermal expansion, dissipative/adiabatic heat production/consumption, 
and here also from the possible dependence of the heat-transfer through the contact surface Fc on the 
delamination itself, and on the possible slot between the bodies if the contact is debonded. 

We consider an elastic response of the adhesive, and then one speaks about adhesive contact (in 
contrast to brittle contact, see Remark 3.5). Within the realm of the literature on (frictionless adhesive) 
contact, in the isothermal case we refer e.g. to [26] in the framework of rate-independent problems. For 
rate- dependent models, we mention [8, 9, 21, 22, 39, 40] (cf. the monograph [48] for further references). 
The anisothermal rate-dependent case has been recently addressed in [10, 11]. The present paper extends 
the analysis in [26] of rate-independent adhesive contact, to encompass inertial, viscous, and thermal 
effects. 

The elastic response in the adhesive will be considered linear, determined by the scalar elastic modulus 
K > 0; cf. Remark 3.4 for a generalization. At a current time, the "volume fraction" of debonded molecular 
links will be "macroscopically" described by the scalar delamination parameter z :Tc [Oi !]■ The state 
z{x) = 1 means that the adhesive is still 100% undestroyed and thus fully effective, while the intermediate 
state < z{x) < 1 means that there are some molecular links which have been broken but the remaining 



Email addresses: riccarda.rossiOing.uiiibs.it (Riccarda Rossi), tomas.roubicek@mff.cuiii.cz (Tomas Roubicck) 



Preprint submitted to Nonlinear Analysis: Theory, Methods & Applications 



April 22, 2010 



ones are effective, and eventually z{x) = means that the surface is already completely debonded at 
x G Fc . As already pointed out in [26] , one needs a specific energy to break the macromolccular structure 
of the adhesive, independently of the rate of this process. Thus, delamination is a rate-independent and 
activated phenomenon, governed by the maximum dissipation principle, and we shall accordingly consider 
a rate-independent flow rule for z. Activating the delamination process in the adhesive contact at a given 
point X S Fc again needs the (phenomenologically prescribed) energy a{x). 

In the thermodynamical context, the energy a{x) needed for delamination is dissipated by the system 
in two ways: one part ai is spent to the chaotic vibration of the atomic lattice of both sides of the 
delaminating surface Fc, which leads "macroscopically" to heat production (cf. also [45, Remark 4.2]), 
while another part og is spent to create a new delaminated surface (or, "microscopically" speaking, to 
break the macromoleculcs of the adhesive). Thus a(x) = ao{x) + ai{x). 

The mathematical difficulties, arising both from the proper thermodynamical coupling and from 
hosting a rate- independent process on Fc , have been already revealed for other inelastic processes in the 
bulk in [44]. The essential ingredient is the satisfaction of the energy balance and, for this, the mentioned 
concept of energetic solutions to rate-independent systems recently developed in [29, 33, 35, 34], and 
adapted to systems with inertia and viscosity in [43] , appears truly essential. 

In Section 2, we set up our model and, in Sect. 3, discuss its thermodynamics and various modifi- 
cations. After making a suitable transformation of the problem using an enthalpy variable instead of 
the temperature, and introducing a suitable weak formulation in Sect. 4, the main existence results are 
presented in Sect. 5, and proved throughout Sections 6-9. For this, in Sect. 6 we set up procedures of 
regularization of the Signorini-type unilateral contact. As we shall observe in Sect. 6, such a regularized 
problem has its own interest. We further approximate it by convexifying some nonlinear terms, and 
setting up a time-discretization procedure in Sect. 7. Hence, we prove fine a-priori estimates. Ultimately, 
a careful passage to the limit is executed in two consecutive steps in Sections 8 and 9. 

2. The model 

Hereafter, we suppose that the elastic body occupies a reference domain 

17 C M"* , d = 2 or 3, bounded and with a Lipschitz boundary dfl. 

We assume that 

with f2+ and il_ disjoint Lipschitz subdomains and Fc their common boundary, which represents a 
prescribed delamination (d— l)-dimensional surface. We denote by v the outward unit normal to dQ, 
and by i^^ the unit normal to Fc, which we consider oriented from il„ to 17+. Moreover, given v G 
VF^'^(17\Fc), (respectively, v~) shall signify the restriction of v to 17+ (to 17_, resp.). We further 
suppose that 

(917 = Fo U Fn , 

with Fo and F^ open subsets in the relative topology of 917, disjoint one from each other and each of 
them with a smooth boundary. 

As state variables, inside 17 we have the displacement u : 17\Fc — > and the absolute temperature 
6 : 17\Fc — > (0, +oo), while on the contact boundary we consider a delamination variable z : Fc — > [0, 1], 
having the meaning of the integrity fraction of the adhesive. Namely, z = 1 (respectively 2 = 0) means 
that the adhesive has full (resp. no) integrity. We denote by 

|u] = u''~|rc ~ Itc — the jump of u across Fc. 

Furthermore, we shall denote by T = T{u, v, 9) the traction stress on some (d— l)-dimensional surface F 
(later, we shall take either F = Fc or F = F^), i.e. 

T{u,v,e) a\^i^ , with a :^ ]D)e{u) + C{e{u)-E9) , (2.1) 

where of course we take as v the unit normal to Fc, if F = Fc- In (2.1), a is the stress (assuming 
Kelvin- Voigt's rheology and thermal expansion, see (3.8) later on). 



2 



To describe various general situations in a unified and simple way, we introduce 
a closed, convex cone K{x) C M'*, possibly depending on x € Fc, 
and assume the boundary conditions on in the complementarity form as 

M h 0, ] 

T{u,u,e) h 0, \ onFc. (2.2) 
Tiu,u,e)iul = J 

In (2.2), y is the ordering induced by the multivalued, cone- valued mapping K : Tc ^ K'', in the sense 
that, for vi,V2 --Tc ^W', 

vi >: V2 if and only if vi{x)~V2{x) € K^x) for a.a. x e Tc- (2-3) 

Likewise, ^ is the dual ordering induced by the negative polar cone to K, in the sense that, for Cij C2 '■ 

Tc K'^, 

Ci ^ C2 if and only if (i{x)-v > C2{x)-v for all v S -ftr(a;), for a.a. x gTc- 



Possible choices for the cone- valued mapping K : Tc ^ K'^ are 

K{x) ^ K = R'^ for a.a. x e Fc, or (2.4a) 

K{x) = {v e M'*; vi'^ix) > 0} for a.a. x E Tc, or (2.4b) 

K{x) ^{v e R''; v-i^^{x) = 0} for a.a. x G Tc- (2.4c) 



In the first case (2.4a), the second of boundary conditions (2.2) translates into T{u,u,9) = on Fc, 
while no constraint on [[uj is imposed. Thus, (2.4a) allows for no interaction of the bodies and 
after a complete delamination. In fact, this model is very simplified because it does not prevent possible 
interpenetration and delamination can be thus triggered, rather unphysically, by mere compression. Nev- 
ertheless, a model like this may be feasible in some situations. In this connection, let us point out that 
the interpenetration after developed cracks is neglected in several crack models used in mathematical 
literature (as e.g. [13, 15, 25]), too. The case (2.4b) yields the standard model of unilateral frictionless 
Signorini contact in the normal displacement at x € Tc- The last case (2.4c) prescribes the normal jump 
of the displacement, variable at a; S Fc, to zero. Thus, it only allows for a tangential slip along Fc. This 
may be a relevant model under high pressure, when no cavity of Fc can be expected anyhow. Such a 
situation occurs, e.g., on lithospheric faults deep under the earth surface. Note that, both in (2.4a) and 
in (2.4c), K{x) is a linear manifold for a.a. x G Tc- As we shall see later, this feature may allow for some 
special benefits. 

Classical formulation of the adhesive contact problem. Beside the force equilibrium, coupled with the 
heat equation inside il\Fc and supplemented with standard boundary conditions, we have two comple- 
mentarity problems on Fc, namely 



QU — div 


{Be{u) + C{e{u)-E0)) = F 




in Q\Y.c, 


(2.5a) 


c^{e)e - 


div {K{e{u), 9)^9) =Be(u):e(u) + 


9CE:e{u) + G 


in Q\^c, 


(2.5b) 


u = Q 






on Sd, 


(2.5c) 


T{u, u, 6) 






on Sn, 


(2.5d) 


{K{e{u),e)^e)v = g 




on E, 


(2.5e) 


[De(u) 4 


■C{e{u)-E9)}iy^ = 




on Scj 


(2.5f) 


H^o 






on Ec, 


(2.5g) 


T{u, u, 6) 






on Ec, 


(2.5h) 


{T[u, u, 6 


>) + Kzluj)-[uj =0 




on Ec, 


(2.5i) 
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i < 






on Ec, 


(2.5j) 


d < fli + ao 






on Ec, 


(2.5k) 


i (d — oq — ai) =0 






on Ec, 


(2.51) 


de A^[oa](z) + [u 


ir 




on Ec, 


(2.5m) 


i(K(e(ri),0)V(?|+, + 


IK(e(7.),0)V0|fJ.j.±+77(H,z)[0 


1=0 


on Ec, 


(2.5n) 


[K(e(w),6l)V6']-i^± = 






on Ec, 


(2.5o) 



where we have used the notation 

Q := (0, T) X f), E := (0, T) x 917, Ec (0, T) x Tc, E^ (0, T) x To, E« := (0, T) x r„, 

T > being a fixed time horizon. In (2.5), F : Q W'- \s the apphed bulk force, wd : Eo ^ M'' the 
prescribed time-dependent boundary displacement, / : E^ — > K"^ the applied traction, while G : Q ^ M. 
and _g : E — > M are some external heat sources. In addition, 

C, D : — > Mfy^^ arc 4th-order positive definite and symmetric tensors, (2.6) 

(i.e. Cijki = Cjife/ = Ckiij, and the same for D), IK = K(e,0) is the positive definite matrix of the heat 

conduction coefficients, and E e M'^^'' is a matrix of thermal-expansion coefhcients. Furthermore, the 

constant k > phenomenologically describes the elastic response of the adhesive. The complementarity 

problem (2.5g)-(2.5i) describes general, possibly unilateral (depending on the choice of the mapping 

K -.Tc ^ M.'^) contact, whereas the adhesive contact results from the complementarity conditions (2.5j)- 

(2.5m). In (2.5k) the coefficient oq (resp. oi) is the phcnomenological specific energy (per area) which 

is stored (resp. dissipated) by disintegrating the adhesive. The overall activation energy to trigger the 

I 1 2 

debonding process in the adhesive is then aq + ai. Note that the term k [uj in (2.5m) is in fact a 
penalization of the delamination condition ~ 0, cf. with the brittle dclamination model (3.18). 

Moreover, A^[o,i] denotes the normal cone to the interval [0,1], i.e. the subdifferential in the sense of 
convex analysis of the indicator function /[o.ij of [0, 1]. Finally, 77 = r]{x, , z) > is a phenomenological 
heat-transfer coefficient^ determining the linear heat convection through Fc. We shall suppose that ry 
depends affinely on the delamination variable z, cf. (5.1e) below. 

3. Thermodynamics of the model and various remarks 

Let us briefly present the thermodynamics of the boundary-value problem (2.5). The underlying 
overall Helmholtz free energy : M'' x M x (0, +00) M has a bulk and a surface part, i.e. 

^iu,z,9)= [ i:^''^^{e{u),9)dx+ [ V'"^([u],z)dS', (3.1a) 

with -0^""^ and ■0*'"'"^, respectively being the bulk and the contact-surface contributions to the specific 
Helmholtz energy. One can identify 

1 0'^ 

0buik^g^^~) ^ -C(e-E6i):(e-E6i) - —M:E-^Po{0) 

= ^Ce:e - 9M:e - ^poie) with B := CE. (3.1b) 

Here ^Ceie is the mechanical part of the internal energy in the bulk, while —tpo{9) is the thermal part of 
the free energy. Hereafter, we shall assume that 

00 : (0, +00) — > R a strictly convex function. (3.1c) 
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The specific contact surface energy , z) is then 

|2 




c I , i„ -I II "'II I anz] dS if [fi;| and 0<2;<1 a.e. on Fn, 

r'''Hv,z)= {Jrc^^ ^ ^ ^ ~~ (3.1d) 



oo otherwise. 



The other underlying ingredient of the model is the overall dissipation rate ^, which also has bulk and 
surface contributions and ^t)""* and ^ namely: 

S(e,z):= / [C(e,i)](dx)- / C^"^e) dx- + / f (z)d5 (3.2) 
Jn Jn\rc Jtc 

where the specific dissipation rate ^(e, z) ~ £^^^^^{e)dx+^^^'-'^ {z)dS is a measure in general, with absolutely 
continuous part determined by the (pseudo)potential of viscous-type dissipative forces in the bulk, and 
a possibly concentrating part, supported on Tq, i.e. 

r'^(e) = 26(e), C2(e) >e:e, f-f(i) = Ci(^) H"^ '[l " ° " 

z I +00 otherwise, 

the latter term representing the potential (and also the specific dissipation rate) of the rate-independent 
delamination process on the contact boundary r^. 

Standardly, one then defines the specific entropy s — s{6, e) by the so-called Gibbs' relation (s, 9) = 
— 'i>g{u,z,d;d) where 'i/'g{u, z,9;9) is the directional derivative of ^ at {u,z,6) in the direction 6. This 
yields the entropy in the bulk as 

^ = -^^(^H,^) = B:e(u) + V^(0). (3.3) 
Further, we shall use the so-called entropy equation 

eh = e''""'(e(w)) - div(j) + G. (3.4) 

Substituting s = W.eiu) — i/'g (0)0, cf. (3.3), into the entropy equation (3.4) yields the heat equation 

c^{9)9 + div(j) = 2C2(e(u)) + eM:e{u) + G (3.5) 

with the heat capacity 

c^{9) - 0<(0). (3.6) 

Hence, postulating the constitutive relation for the heat flux 

J :=-K{e{u),9)V9, (3.7) 

i.e. Fourier's law in an anisotropic medium, one obtains to the heat equation in the form (2.5b). 

Similar, but simpler thermodynamics can be seen also on the contact boundary by involving ip^^'-'^ 
and , As (3. Id) is independent of temperature, the "boundary entropy" =—^V^"'^^ is simply zero, 
and the corresponding entropy equation reduces to = £^^™^{z) — [j] (as an analog of (3.4)), which then 
results in (2.5o). Incorporating the analog of the phenomenological law (3.7), we arrive at (2.5n). 
Momentum equation. As in Kelvin- Voigt rheology, the total stress a is postulated as 



-=^C2(e(^)) + ^^^"'^(e(^),e), (3.8) 



which just gives a from (2.1). From Hamilton's principle, generalized for dissipative systems as in [6] and 
with the specific kinetic energy ipj-itp, one then obtains the equilibrium equation (2.5a). For later use, 
we introduce the indicator functional Ik associated with K -.Tc ^ K'', defined on L^(rc;M'') by 

Ik{v)^ [ lK^^){v{x))dS{x) for all u e L2(rc;M''). (3.9) 
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Wc point out that the complementarity conditions (2.5g)-(2.5i) may be reformulated as the subdifferential 
inclusion 

dlKiluj) + T{u,u,0) + Kzluj 3 inSc, (3.10) 

featuring the (convex analysis) subdifferential OIk ■ i^(rc;M'*) =^ L^(rc;M'^) of the indicator functional 
Ik introduced in (3.9). 

Evolution of the delamination parameter (a flow rule). Finally, we consider the following differ- 
ential inclusion for the inelastic evolution of the parameter z 

9Ci (2) +5.i^™'''(u,z) 9 inSc, 

which in the adhesive case results in 

dl(-oo,o]{z) + dI[o^i]{z) + ^k\IuI\'^ - ao - ai 3 in Ec- (3.11) 

It is immediate to check that (3.11) is a reformulation of (2.5j)-(2.5m). 

The entropy equation (3.4) is designed to balance the total energy, i.e. the sum of the kinetic energy 
integrated over i}\Tc with the overall dissipated energy (i.e., S from (3.2) integrated in time), and with 
the bulk internal energy 

e'""""(e, 9) := V'''""' + ds ^ ^Ce{u):e{u) - Vo(«') - OM:E + 0{M:E + i''„{e)) 

= h{9) + ^Ce{u):e(u) with h{e) := ei''„{e) - i^oi^); (3.12) 

we convene to refer to h as the enthalpy, see also [44, Sect. 2]. One can then derive the total energy 
balance: 



dt 



^ q\u\^ + -<Ce{u):e{u) + h{e)Ax 
n\Tc 2 

kinetic, clastic, and thermal energies 
|2 



aozdS = / G + F-uAx + [ gdS+ [ f-iidS . (3.13) 
Jn Jr JFn 



mechanical energy power of bulk heat power of surface heat 

in the adhesive and mechanical load and mechanical load 

Assuming Oq > 0, G > a.e. in ft, and <? > a.e. in d^l, we can rely on the fact that 9 > a.e. in 
(proved later in Theorem 5.1) and, using (3.4), we derive the Clausius-Duhem inequality: 

Remark 3.1 (Partly linearized ansatz). An important feature is that, as a consequence of the partly 
linearized ansatz (3.1b), the mechanical and thermal variables are additively separated in (3.3), which 
makes Cv in (3.5) independent of u, and thus makes mathematical analysis easier. 

Remark 3.2 (Non-homogeneous boundary conditions). We could supplement (2.5a) with non-homoge- 
neous, Dirichlct boundary conditions on To, i.e. impose 

u = uiu on Sd, (3.15) 

for some prescribed time-dependent loading ujd : [0,T] — s- iJ^/^(rn). The analysis we are going to 
perform in the case of homogeneous Dirichlct conditions can be carried over to the case of (3.15) by 
arguing as in [45], and thus recurring to the additive split u{t) ~ u{t) + UD(i) for almost all t G (0,T), 
with u : [0,r] Wr;^(fi\rc;M'^) and uu : [0,r] W^''^{n\Tc;M.'^) an extension of ujb to fl. It was 
observed in [45] that, if Fc H Fd = 0, one can assume that [['«D(i)l = 0, whence lu{t)} = lu{t)} for almost 
all t G (0, T). This allows for a reformulation of the problem in terms of the unknown it, hence reducing 
the analysis to the case with homogeneous Dirichlct conditions. 



6 



Remark 3.3 (Heat-transfer contact conditions). Note that the transient conditions (2.5n)-(2.5o) on Tq 
for the heat equation can equivalently be written as two Robin-type conditions 

K(e(u),0)V0|+^-;.±+77(H,^)0|+^ =77(H,^)e|F, -^aii on Sc, (3.16a) 
K{e{u),0)Ve\^^;y^+v{luj,z)e\^^=ij{luj,z)e\+^-^a,z on Sc, (3.16b) 

where we have highhghted the unit normals from i7_ to 51+ and from to n_. This reveals that 
the heat generated by delamination is distributed with proportions ^ and ^ into the two subdomains 
adjacent to Tq. In principle, we could also consider a contribution from the Stefan-Boltzmann radiation, 
which would then result in the condition 

i(K(e(u),0)V0|+^+]K(e(u),0)V(?|rJ-t^±+%(H,2)[0] + 7?i( [h] , z) [^^j ^^^.^^ (3^^^) 

with ?/o([w],-)i '7i(II'"1j') ^ affinc. In fact, this would lead to a modification of the present analysis 
which is quite routine. Thus, we shall not scrutinize this generalization here. Let us also remark that, in 
alternative to the dependence of 77 on |m] , a dependence on the normal stress is sometimes considered, 
cf. [1]. However, it seems difficult to adapt the present multidimensional analysis to that case. 

Remark 3.4 (Elastic response in adhesive). One can easily imagine a positive-definite dxd-matrix in 
place of K, which would more properly describe the phenomcnological clastic response of the adhesive. 
The related analysis would be just a standard modification of the presented one. 

Remark 3.5 (Grifiith concept). The classical concept of delamination is based on the Griffith criterion 
[23], phenomenologically prescribing the amount of energy a (in J/m^, in 3-dimensional situations) needed 
to delaminate the surface, independently of the rate of the process. The classical Griffith-type approach 
considers the adhesive inelastic and one speaks of a brittle delamination. Our adhesive contact problem 
can be viewed as a regularization of this brittle delamination, and in fact makes mathematical analysis 
and numerical implementation easier. It has its own interpretation and many applications, and it is thus 
often considered as the original problem, cf. [8, 9, 10, 11, 26, 40]. In the quasistatic isothermal case, it 
has been proved in [45] that the adhesive contact approximates the brittle delamination as the elastic 
modulus K — >■ 00 in the framework of the so-called energetic solution concept. Furthermore, in [32], any 
energetic solution to the brittle delamination has been proved to be of "Griffith-type" in the sense that 
z indeed takes either the value 1 or the value 0. 

Remark 3.6 (Engineering models). In the engineering literature, the GriHith-type delamination on 
a prescribed so-called "weak surface" is a quite accepted concept (for example in the framework of 
the so-called Finite Fracture Mechanics), although it is often combined with the heuristically devised 
stress criterion, which in some situations seems to provide a better understanding of the initiation of 
the delamination process, cf. e.g. [27, 28]. The initiation of the delamination process can sometimes be 
triggered by another crack approaching the weak surface, according to the classical so-called Cook-Gordon 
mechanism [14], which has been confirmed experimentally. The present form of the activation energy 
a{x) may typically correspond to crack growth in a pure fracture mode (e.g. Mode I). Nonetheless, it is 
believed that the present approach can be extended to a generalization of this form, in order to cover 
more complex phenomcnological engineering models, working with the so-called "fracture mode mixity" 
which reflects the character of the load (the ratio of its shear and normal components) on the crack tip. 

Remark 3.7 (Constant heat capacity). The special case ipo{0) = co9hi{9/6Q), with cq > and 60 > 
constant, would give Cv{0) = cq in (3.6). However, this case is not within the scope of our analysis, since 
Cv(-) does not have a compatible growth, cf. (5.1b). 

Remark 3.8 (Brittle delamination model). Let us now briefly comment on the model for brittle delam- 
ination with thermal effects which would result from the above derivation. As in the case of adhesive 
contact, we focus on its classical formulation, which couples the momentum equilibrium equation (2.5a), 
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the heat equation (2.5b), the boundary conditions (2.5c)-(2.5f) and (2.5n)-(2.5o) with the two fohowing 
complementarity problems on Tc- 



M ^0 


on Ec, 


(3.18a) 


T(w, u,e)hO wherever z(-) = 


on Ec, 


(3.18b) 


T{u,u,0)iuj =0 


on Ec, 


(3.18c) 


z[u}=0 


on Ec, 


(3.18d) 


z <0 


on Ec, 


(3.18e) 


d < ai + ao 


on Ec, 


(3.181) 


z{d — flo — ai) = 


on Ec, 


(3.18g) 


rfe%,i](^) + a.J(H,z) 


on Ec. 


(3.18h) 



Indeed, (3.18a)-(3.18c) and (3.18d)-(3.18h) respectively correspond to (possibly unilateral) contact and 
activated delamination. Note that the penalization terms in (2.5h)-(2.5i) and k/2||m]P in (2.5m) 
are no longer present, and the delamination constraint (3.18d) is enforced by the second subdilFerential 
operator in (3.18h), featuring the indicator lunction 

J(^,^)=/{..=o}, i.e. •^K^)={+^ iJtherwis; ^^'^^^ 

In lact, as function of the two variables v and z J is nonconvex, but separately convex. Hence, the 
subdifferentials of the convex functions J(., z) and J(v, ■) are well-defined, and in particular 9z J([[u]], z) 
is given by 

r if 2 7^ and lu} ^ 0, 
9;, J([u] , 2) = < if z = and M 7^ 0, 
[ R if z 7^ . 

As we have already mentioned, existence for the (global) energetic formulation of the brittle delam- 
ination problem in the isothermal quasistatic case has been proved in [45]. In contrast, the analysis of 
the corresponding thermomcchanical model given by (2.5a-f,n-o)-(3.18) is for the moment being an open 
problem. The main difficulties attached to this problem are related to the presence of two multivalued 
operators in (3.18h), and in particular to the essentially nonconvex character of the nonlinearity (3.19). 

However, taking into account (3.18h), we clearly identify a drawback of the differential formula- 
tion (3.18d-h) of brittle delamination. Indeed, in this framework any driving tendency towards delami- 
nation is smeared out if < z < 1, because then the driving force is c? = < ao+ai, and necessarily, by 
(3.18g), we have z = 0. The adhesive contact problem shows a similar behaviour if k — >■ cxj. 

4. Enthalpy transformation and energetic solution 

Throughout the paper, we shall adopt the notation 

W"r;f (n\rc;K'^) ^ {v e W^'^{n\rc;M.'^) : w = on To} , 
Wr;^^(f^\rc;R'^) = {w e H/i^^^^^^p^.jgd) . „ = q on Tc} . 

Furthermore, in the case K{x) is a linear subspace of R'^ for almost all x £ Fc, we shall use the notation 

W^^{n\ra;M.'^) ^{v£ W^'^{n\rc;M.'^) ■. lv{x)j e K{x) for a.a.x- e Fc}. 
We shall also extensively exploit that, for d < 3, 

W^-^{n) C LP{n) continuously for 1 < p < 6, 

u^u\r: W'^'in) ^ H'/\T) C L'"(F) / ^^^^tinuously for 1 < m < 4, (4.1) 

I compactly for 1 < m < 4, 
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with T = dn, OT T ^ Tc, or r — r^. The same embeddings hold for the Sobolcv space Vl^^''^(r^; M^) 
of vector- valued functions. Finally, we shall denote by (•,•) the duality pairing between the spaces 

Wp;f(17\rc;M'*)* and VKp;^^(rj\rc; M'*). 

The analysis of the nonlinear heat equation (2.5b), featuring the quadratic coupling terms with the 
momentum balance equation (2.5a), calls for rather sophisticated techniques and suitable working as- 
sumptions. In particular, one may impose some conditions either on the growth of K(e,-) (cf., e.g., 
[19] for the analysis of a similar nonlinear heat equation in some phase transition model), or on the 
growth of Cv (cf., e.g., [42, 44] and, more specifically, [17, Sect. 5. 4. 2] for contact problems in thermo- 
viscoelasticity) . Under the latter kind of assumptions, the Galerkin approximation method for proving 
existence of solutions could serve quite effectively, cf. [42] . 

On the other hand, system (2.5) hosts the delamination rate-independent process on Tc- Hence, the 
Rothe method (i.e. the implicit discretization in time) seems more natural for the analysis, see e.g. [35, 
29, 24]. In turn, the nonlinearity Cv(-) makes it technically difficult to implement such a discretization 
method. This problem can be circumvented by rewriting the original PDE system (2.5) in terms of the 
enthalpy, instead of the temperature, as e.g. in [42]. 

Namely, we introduce the so-called enthalpy transformation, setting 



^^h„ie) := / Cv(r)dr. 
Jo 



(4.2) 



Thus, ho is a primitive function of Cv, normalized in such a way that /io(0) = 0. In view of (3.6) and (3.12), 
we have 



h'i9) = (eyj'.ie) - MO))' = + V^(^) - i'M = ^V^'(^) = cM = Kio), 



(4.3) 



hence ho differs from h just by a constant, namely "00(0). Furthermore, thanks to (3.1c), Cv is strictly 
positive and hence ho is strictly increasing. Thus, we are entitled to define 



e(^) 



_ / ho\d) if i?>0, 
] if < 0, 



K(e,e(t^)) 



(4.4) 



where Hq here denotes the inverse function to h. 

Taking into account (4.4), as well as the subdijferential reformulations (3.10) and (3.11) of the com- 
plementarity problems (2.5g)~(2.5i) and (2.5j)-(2.5m), respectively, the PDE system (2.5) turns into 



gu-div{]n)e{u)+Ce{u)-Me{^)) =F 1 
^ - div(3<:(e(w), d)\7d) = ID)e(u):e(u) + e{d)S:e{u) +G J 
u = 

{X{e{u),d)Vd)i' = / 
7{u, ii, 1?) = g 

lBeiu)+Ce{u)-Me{^)P^ = 

diKiM) + T{u,u,e) + Kzluj 3 

91{-oo,o]iz) + dI[o^i]{z) + ^k\IuI - flo - ai 9 

i {X{e{u),d)^d\+^ + X{e{u), z^)Vi?|f J -z/i + 77(M , z)[e(z?)l = 

|3C(e(u),i9)Vi?]]-i/± = -aiz 



in Q\J^c, 

on So, 
on Sn, 
on E, 



on Sr 



(4.5a) 

(4.5b) 
(4.5c) 
(4.5d) 



(4.5e) 



where 



7{u,v,^) := T(w,u,e(z?)) ^ [Be{v) + Ce{u) - 
where again we take as i' the unit normal to Fc j^*, if F = Fc- 



(4.5f) 
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Data qualification. Hereafter, the problem data F, G, J , and g shall be qualified by 

F e Li(0,T;L2(rj;R'^)); (4.6a) 

ri¥i'i(o,r;i4/3(r„;R3)) ifd = 3, 

|l¥i^i(0,r;ii+<^(rN;M2)) for some e > 0, if d = 2; 
GeLi(g), G>Oa.e. inQ; (4.6c) 
g e L^(I]), g > a.e. in S . (4.6d) 

The energetic formulation associated with system (4.5) hinges on the following energy functional $ 
(which is in fact the mechanical part of the internal energy (3.1a)), and on the dissipation potential 3? 



<S>{u,z):= f ^Ce{u):eiu)dx + lK{luj)+ f (|z| [u] |' + /[o,i] (z) - aoz) d5 (4.7) 



lce{u):eiu)dx + lK{luj)+ [ ' 
n\rc ^ JTc 

a\\z—z\dS if z < z a.e. in Fc, 
oo otherwise. 




For notational convenience, we also set for all v e i (51) 
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We are now in the position of introducing the notion of weak solution to system (4.5) which shall be 
analyzed throughout this paper. The reader is referred to [44, Prop. 3.2] for some justification of the 
energetic solution concept in the framework of general thermomechanical rate-independent processes, in 
particular for the proof of the fact that energetic solutions are also conventional weak solutions whenever 
z is absolutely continuous. 

Definition 4.1 (Energetic solution of the adhesive contact problem). Given a quadruple of initial data 
(wq, uo, ^Oi ^o) satisfying suitable conditions (cf. (5.4) later on), we call a triple (u, z,i?) an energetic 
solution to the Cauchy problem for (the enthalpy reformulation of) system (4.5) if 

u G W''^{0,T;W^i^{n\Tc;R'')), (4.9a) 
u&W^^°°{0,T;L'^{n;R'^)) ii g > 0, (4.9b) 
z e i°°(I]c) n BV([0,T];ii(Fc)) , (4.9c) 

BV([0,T];M^i''^ (0\Fc)*) J ^r every i s r < , ^4.jaj 

with r' denoting the conjugate exponent ^^r^ of r, and the triple (w, z, ■&) complies with: 

(i) (weak formulation of the) momentum inclusion, i.e.: 

{uj h on Sc, and (4.10) 
gu{T)-{v{T)-u{T)) dx + / (De( u) + Ce(u) — B0(?9)) :e{v—u) — gu-(v—uj dxdt 

Jq 

+ [ Kz[uj-lv-ujdSdt> [ guo-{v{0)-u{0))dx+ [ Fiv-u) dxdt + [ /•(w-u)dS'di (4.11) 

for all V in L^{0,T;W^'^{n\Tc;l^'^)) with fvj ^ on Sc and, if g > 0, also in W^'^{0,T; L'^{n;R'^)), 

(ii) total energy inequality 



T^,MT))+<P{uiT),z{T))+ lyiT)dx<T^.Juo)+<f{uo,zo)+ jj^dx 

Gdxdt+ / gdSdt (4.12) 



/ F-udxdt + 


I f-udSdt+ [ 


Iq 
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(iii) semistability for a. a. t G (0,T) 

V5€L°°(rc): ^{u{t),z{t)) <^{u{t),S) z{t)) (4.13) 

(iv) (weak formulation of the) enthalpy equation: 

d{T)w{T) dx+ [ X{e{u), i})\/^-Vw - dw dxdt + [ vi H , [6(1?)] {wj dSdt 

n JQ JT,c 

(©e(w):e(w) + e(z9)B:e(w))wdxdt+ / ^'Irc+^lrc ^surf^^^.^^^ 

G'wd.Tdt+ I gwdSdt+ [ dow{0)dx (4.14) 



Jn 

for all w e C°([0,T]; M^i''''(17\rc)) n Vl^i''''(0, T; L'-'(f7)), where i?o := ho{9o) and is a measure 
(=heat produced by rate-independent dissipation) defined by prescribing its values for every closed 
set of the type A := [ti,t2]xC C [0,T] x T\r as 

f, , I ai\z{ti,x)-z{t2,x)\dS if z(-,x) nonincreasing on [ti,t2] for a.a. zeC, -,r^ 
[ +00 elsewhere, 

(v) and the remaining initial conditions (in addition to ?i(0) = wqi already involved in (4.11)), i.e. 

u(0) = Mo a.e. in il, z(0) = zq a.e. in Fc, i9(0) = a.e. in VI. (4.16) 

Remark 4.2. The inequality (4.12) is the integrated (inequality) version of the total energy bal- 
ance (3.13). It is immediate to check that, for every closed set of the type A := [ti,t2\y~C C [0,T] x Fc, 
^™'^(A) coincides with Varg^(z|c; [^1,^2])- Now, relying, e.g., on [20, Prop. 1.3.10,Thm. 1.5.6], one can 

verify that formula (4.15) indeed defines a non-negative Radon measure on Ec. Subtracting (4.14) tested 
by 1 from (4.12) reveals the mechanical energy inequality: 



T^r.{'^{T))+'^{u{T),z{T))+ / Be(M):e(u)dxdi + Var3j(z;[0,T]) 

JQ 

<T^.^{uo) +^{uo,zo) + [ F-u~e{d)M:e{u)dxdt+ f f-udSdt. (4.17) 

JQ "'Sn 

In particular, z{-,x) must be nonincreasing on [0,T] for a. a. x € Fc, otherwise Var3j(z; [0,T]) = 00 and 
(4.17) cannot hold. 

5. Main results 

Assumptions. Hereafter, we shall denote by the symbols C, C most of the (positive) constants occurring 
in calculations and estimates. We suppose that 

Cv : [0, +00) — > continuous, (5.1a) 

> w > J^, ci > Co > V6I e M+ : co(l+6l)""i < c^(6i) < ci(l+6i)"i-\ (5.1b) 

K : R'''"^ X M ^ M''^'^ is bounded, continuous, and (5.1c) 

inf X{e,d)C.^ = k>0, (5. Id) 

and that ri{x, v, •) is a non-negative affine function of the delamination parameter z e [0, 1], i.e. 

ri{x,v, z) ~ i'ii{x,v)z + rio{x,v) for 771 , 770 : Fc x M'^ — > Caratheodory s.t. 

3C^>0: V(a;,t.) eFcXM'^ |77o(a;,7;)| + |7/i(x,^;)| <C;^(|t.|4/3 + l); ^^'^""^ 
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in fact, the above growth condition for the functions riQ{x,-) and 7yi(x,-) is not optimal and could be 

slightly improved, as one can deduce from the proof of Theorem 6.1 in Sect. 8 later on. It is immediate 
to deduce from (5.1b) that 

3C^, Cl > Vw e M+ : CHw^'"^^ - 1) < Q{w) < Cl(y^l'^ - 1) . (5.2) 

Moreover, it follows from (5.1c), (5.2), and the definition (4.4) of 3C that 

3Cgc >0 Ve, CeM'^ : |3C(e,^)eCI <Cgc|?||C|. (5.3) 

Finally, we impose the following on the initial data 

uo e Wl^li^XV^-^'') , M ^ on Sc, (5.4a) 

uo e i^(^^;M'') if £0 0, (5.4b) 

zo G i°°(rc), < zo < 1 a.e. onTc , (5.4c) 

6loei"(r2), 6*0 > a.e. infi. (5.4d) 

Theorem 5.1 (Existence for the adhesive contact problem). Let us assume (4.6), (5.1), (5.4) and 

(i) if g = (such a case is sometimes referred to as quasistatic let also 



[W^'^{0,T;L^+''{n;R^)) for some e > 0, if d ^ 2; ^ ' 

M"^'^ {d^^r\Y^)>Q, ''-I (9f}_ n To) > 0, (5.5b) 

with denoting the (d—l)- dimensional Hausdorff measure, or 

(ii) if 9 > l^t o,lso 

Q > {) and K{x) is a linear subspace ofW^ for a.a.x e Tq. (5.6) 

Then, there exists an energetic solution (u, z, -d) to the adhesive contact problem with the additional 
regularity 

u€W^''^{0,T-w]f{n\Vc-R'^)*) ifg>0. (5.7) 
Furthermore, in both cases g > and g = 0, the positivity of the initial temperature 

inf 9q =:9*>0 (5.8) 

implies inf .jjgQ = inf(t ,j,)gQ 0('!9(<, x)) > 0; in particular, 9 is a.e. positive on Q. 

Theorem 5.1 shall be proved in Section 9 by passing to the limit in some regularized problem (where 
the contact conditions on Tc are penalized), which we shall present in Sect. 6. In turn, existence for the 
latter problem shall be proved in Section 8 by passing to the limit in a further approximation scheme, 
constructed in Section 7 by a regularized semi-implicit time discretization. 

Remark 5.2. In Theorem 5.1 we distinguish the cases g > and g = 0, because in the latter case we are 
able to prove existence for a far larger class of cones yielding the unilateral constraint on the displacement, 
in particular the Signorini conditions. This stems from the fact that the analysis of the momentum 
equilibrium equation in which inertia interacts with Signorini boundary conditions is remarkably difficult. 
It has indeed been an open problem for a long time and only very recently, in [37, 38], existence results 
have emerged for the dynamical viscoelastic equation with Signorini contact conditions, in the one- and 
three-dimensional case on unbounded domains. Such results have been proved with very sophisticated 
Fourier analysis techniques. In the one-dimensional framework of [38], it has also been obtained that the 
solutions comply with the energy balance. 
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Remark 5.3. Under (5.6), the qualification v G W^^'^{0,T; L'^{n;M.'^)) for the test functions in (4.11) 
might be relaxed to 

we W^'^{0,T;Wj/{n\Ta;R'^)*). (5.9) 
Indeed, thanks to (4.10) and to the linearity of K{x) for almost all x G T^, the function u fulfill- 
ing (4.11) is such that u e L^{0,T;W]^^{n\Tc;l&'^)). Now, the spaces L^{0,T;W^^{n\Tc;'^'^)*) and 
L^(0, T; Wj^'^(r2\rc; K'')) arc in duality. Hence, (5.9) is sufficient to give meaning to the term J^u-vdx. 
A similar extension will also apply for the test functions of (6.7) below. 



6. Regularization 

We shall approximate (the enthalpy-reformulation) of the adhesive contact system (4.5) by penalizing 
the contact condition Ju] ^ 0. This is a well-established routine in the analysis of contact problems, see 
e.g. [17]. We should emphasize that the penalized problems themselves have their own practical usage 
because they allow, first, for combination of inertia and the unilateral- type elastic contact condition 
and, second, for a more physical interpretation of the coupling through the heat-transfer coefficient, 
cf. Remark 6.2 below. 

Thus, we shall replace the subdifferential operator OIk in the differential inclusion (3.10) (equivalent 
to the complementarity problem (2.5g)-(2.5i) on Sc), with its e-Yosida regularization (see, e.g., [2, 12, 3]). 
We recall that the e-Yosida approximation of the indicator functional Ik is the lower semicontinuous, 
convex, and Frcchet diffcrentiable functional given by 

/|, :L^(rc;R'^)^ [0,+oo) given by If^iv) ^ ^min\\v - w\\l,^^^.^^,y, (6.1) 

cf. definition (2.3) for the ordering ^. We point out that Mosco-convcrges to Ik in the L^(rc;R''); 
see, e.g., [2, § 3.3] for the definition of Mosco-convergence and [2, Thm. 3.66] for the link with Yosida 
regularizations. In particular, this entails that 

Ve^v in L^iVc^R"^) liminf /^-(we) > //f(w). (6.2) 

The e-Yosida regularization of OIk is the Frechet derivative (I^)' : L2(rc;M'^) ^ ^^(rc;^'^) of the 
functional It is well known that 

(/|,)' = i(Id-PK), (6.3) 

where Id : L2(r^;Rd) L^^p^.j^d^ jg identity operator and : L'^{Tc;'M.'^) L'^iTc,'!^'^) is the 
projection associated with the multivalued mapping K : Tc ^ K''. For later use, we recall that, being 
Pk a contraction on L^(rc;M''), there holds 

ll(^K)'(^^)llL^(rc;K'') < ^IK'llL^(rc;K'') ^r all V e L\ra;R''). (6.4) 

Hence, we shall consider the following regularized conditions on Fc, where (3.10) is approximated by 
Yosida regularization: 

|ID)e(u) + Ce(u)-Be(i9)lzy± = 
Kzlul + {II,y{lul) + T{u,u,e)^0 

dI{-oofi] (z) + 9/[o,i] (z) -f iftl |u] I" - flo - ai 9 
^{X{e{u),^)\/^+^ + X{e{u),d)Vd\^^);y^+7j{lu},z)iem=0 

The resulting regularized stored energy is then 

<i>e{u,z):^ [ ^Ceiuy.e{u)dx + II,{luj)+ [ (^z\[ul\' + I[o,^]{z) - aoz) dS. (6.6a) 
Jn\rc ^ Jrc^^ ' 

The main result of this section ensures the existence of energetic solutions to the initial-boundary value 
problem for the adhesive contact model supplemented with the regularized contact conditions (6.5). 



on Ec. (6.5) 
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Theorem 6.1 (Existence of energetic solutions to the regularized problem). Under assumptions (4.6), 
(5.5), (5.1) and (5.4), for every £ > there exists a triple [u^, z^^-d^) as in (4.9), and such that, in 
addition, 

u, e w^^\Q,T-wl-^{n\T^-m.''y) ifg>o, 

which solves the energetic formulation of the Cauchy problem for system (4.5a-d) and (6.5), namely the 
initial conditions (4.16) hold, as well as 

(i) the (weak formulation of the) momentum equation: 

u,)+<Ce{u,)^MQ{-d,)):e{v)- Qu,-idxdt+ [ (kz4u,] +(/^)'( [u,] )) • [t;] d5di 

+ / gu,{T)-v{T)dx ^ I guQ-v{0)dx+ [ F-vdxdt+ [ f-vdSdt (6.7) 
Jn Jn Jq JEn 

for allv in L^{0,T;W^'^{n\rc;^'^)) and, m the case g>0, in W^'^{0,T;L^{n;R'^)), 

(ii) the total energy inequality (4.12) with and {u^, z^,'d^) in place of ^ and {u,z,{}), 

(Hi) Zg, complies with the semistability condition (4.13) for a. a. t £ (0,T), again with $5 and (u^j^e), in 

place of $ and (it, z) 
(iv) the weak formulation (4.14) of the enthalpy equation. 
Furthermore, if (5.8) holds, then 

inf ^Jt,x)>0. (6.8) 

e>0, {t,x)eQ 

Remark 6.2 (Signorini-contact case). Let us point out that, if A' : Fc =t M'' is of the type (2.4b), 
i.e. corresponding to unilateral frictionlcss Signorini contact on Fc, the e-Yosida regularization of OIk 
is given by (/^)'(ue) = — ^[[weln '^ith [[we]ln = [wel'J^^ and with (•)^ = — niin{0, •}, and the second of 
(6.5) reduces to 

Kz[u,j -^[u,j;^iy^ +Tiu„u„9,) =0 on Ec. 

It follows from the above relation that, for fixed e > 0, in the case and z{t,x) > one can express |{u| 
as a function of the traction stress T{u,u,6). This also holds for in the case for z{t,x) = and 
Tn{u,u,9) > 0, while the tangential stress Tt{u,u,9) = because there is no friction. Let us suppose 
that the heat-transfer coefficient ri{-,z) vanishes if there is no contact, i.e. on the set fuja < 0}: 

then, by continuity, ri{-,z) vanishes also on {[uj; \u\n < 0}, and hence by substitution one can express 
the heat-transfer coefficient as a function of the normal stress and of z. We point out that the mentioned 
condition on ?7(-, z) is, to some extent, a natural assumption, also advocated in the engineering literature, 
cf. e.g. [46]. Such an approach does not seem mathematically amenable for the multidimensional non- 
penalized Signorini problem; for c? = 1 we refer to [1]. 

Scheme of the proof. The proof of Theorem 6.1 shall be developed in the next subsections by pursuing 
the following scenario. First, in Section 7 we shall devise a semi-implicit time discretization (with a 
further regularization in the momentum equation), and prove existence of solutions to the time-discrete 
problem. Next, in we shall derive refined a priori estimates, enabling us to perform the limit passage as 
the discretization time-step r goes to in Sect. 8. In this way, we shall conclude the proof of Theorem 6.1. 



7. Semi-impHcit time discretization 

We perform a semi-implicit time-discretization using an equidistant partition of [0,T], with time-step 
T > and nodes := fcr, fc = 0, . . . , A',-. 

Hereafter, given any sequence {</>^}j>i, we shall use the following notation for the backward difference 
operator and its iteration by, respectively, 

T 
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Secondly, we recall the notion of picccwise constant and piecewise linear intcrpolants: for a given AV-tuple 
{^rlfcli C 25, (23, II • lis) being some Banach space, the left-continuous piecewise constant interpolant 
br : (0,T) — > !B, the right-continuous piecewise constant interpolant &^ : (0,T) — > 25, and the picccwise 
linear interpolant : (0,T) 23 of the elements are the functions respectively defined by 

Mt)^b';, b,{t)=b^-\ br{t)^^—^bl+^-^^bl-' forte (7.2) 

T T 

Thirdly, we shall denote by tr and by the left-continuous and right-continuous piecewise constant 
interpolants associated with the partition, i.e. \.r{t) = if t^~^ <t <t^ and t^{t) = t^~^ if t^~^ <t<t^. 
For later use, we recall the following elementary inequalities for all t G [0,r] 

\\br{t)U < \\MmS + \\bMB=\\Mt)h + \\brit-Ut))h, (7.3) 

\\br{t) - bAt)\\^ < r Wh'MUds. (7.4) 
Time-discrete problem. Wc approximate the data F, / by local means, i.e. setting for all fc = 1, . . . , 

- r F{s) ds , f'^ - f{s) As , 

and consider the interpolants Fr, /,-, and fr of the /v-tuples {i^r l^lii {/tIa^i- view of (4.6a)-(4.6b), 
the following estimates and strong convergences hold as r — > 

^ \ inL-P{Q,T-L'^/\Q.-W^))ioT B\\l<p<oo if g = 0; ^ ' 

3C>0 Vt>0 : ||/rllL=(0,T;L-i/3(rN;R'')) - '^ll/ll-L°=(0,T;L*/3(rN;R'')) ' 

7,^/ in iP(0,r;i'*/^(rf,;R'')) for aU 1 <p < oo as 0, (7.5b) 

3C > Vr > : ||/T||Li(0,T;L'i/3(rN;R'')) - 2||/||Li(0,T;L*i/3(rN;R'')) • 

Furthermore, we shall approximate G and g with suitably constructed discrete data {G^}^^^, {g'^}^Zi 
with 

e M^^'2(17)*, g'!^ e H^/^{dny for ah = 1, . . . , i^^, (7.6a) 

and such that 

Gr^G in Li(Q), ^ in L^i^) as t ^ , (7.6b) 

and approximate the initial datum uq with a sequence {uo^t} C Wj!;^''(ri\rc; M'') (with 7 > max{4, ;jrr} 
as assumed in Problem 7.1) such that 

lim /F||e(uo,r)||L-'(n.Rd) = 0, wo,r ^ in W^-^{n;W^) as r ^ 0. (77) 

We are now in the position of formulating the time-discrete problem, which wc again write in the classical 
formulation for notational simplicity. 

Problem 7.1. Let 7 > max{4, J^}. Given 

find {(w^^, 1?^^, fulfilling for k = 1,...,Kt the equations in Q\Tc 

gB^tU^^, - div(]D)e(Dtii^,) + Ce(w^,)-Be(i9^,) + T|e(ii^,)|^"'e(w^,)) = (7.9a) 
DtT?^^, - div(2C(7?^,, e(u^;))V79^;) - ^:^Pe(Dt7.^^,):e(Dtu^;) + e(i9^,)B:e(Dt^i^,) + G^, (7.9b) 
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with the boundary conditions 

on To , (7.10a) 

onr^, (7.10b) 

on on, (7.10c) 

on Tc , (7.11a) 

on Tc , (7.11b) 



on Tc , (7.11c) 

onTc, (7.11d) 
on Tc ■ (7. lie) 

Remark 7.2 (Semi-implicit discretization). The value u^^^ the level /c — 1 in (7. lid) makes the above 
scheme semi-implicit, not just fully implicit as it would be in the case u^^ were in place of u^^T^ in (7. lid). 
This makes the proof of Lemma 7.4 easier. 

Remark 7.3 (Regularization). Like in [44], a regularizing term T|e(u)|'''~^e(M) was added to the mo- 
mentum equation in the bulk and to the corresponding boundary/contact conditions, too. Its role is 
to compensate the growth of the right-hand side terms in the momentum equation, cf. the proof of 
Lemma 7.4. Moreover, with the aim of obtaining some suitable semiconvexity property of the approxi- 
mate stored energy (cf. Lemma 7.5 below), we have also introduced the monotone terms 

r"z and T''(l-h|[u]|^)^"^[w], with 4 < /i < 5 and a, /3 e (0, 1) (7.12) 

in the differential inclusion for the delamination parameter and in the boundary conditions for w on Fc, 
respectively; see Lemma 7.5 for some further specification of the exponents a and /?. 

Lemma 7.4 (Existence of weak solutions to Problem 7.1). Under the assumptions of Theorem 6.1, for 
every fc = 1, Kr there exists a triple (w^^, z^^, i?^V) € iyr;7(rj\rc; K'') x L°°(rc) x W^^^{n\Tc), fulfilling 
the weak formulation of the boundary value problem (7.9)-(7.11). Moreover, > a.e. in il. If, in 
addition, (5.8) holds, then there exists some constant x* > (cf (7.23) below) such that, for sufficiently 
small T, 

'!^er>X*>0 a.e. inn for every k = I,..., Kr. (7.13) 

Proof. The existence of a weak solution to the boundary value problem (7.9)-(7.11) can be proved relying 
on the standard theory of pseudomonotone set- valued operators (see e.g. [41, Chap. 2]). In particular, we 
may apply Leray-Lions type theorems. Indeed, the strict monotonicity of the main part of the operator 
which comes into play in the weak formulation of problem (7.9)-(7.11) derives from the presence of the 
term r|e(u)|'''~^e(u). The latter counteracts the quadratic nonlinearity in e(u) of the dissipative heat 
source in (7.9b). 

We now show that this operator is coercive w.r.t. the norm ofWp'^ {n,\Tc; K'^)xi°°(rc)x VF^'^(f2\rc)- 



(De(DtM^,) -f Ce{u^^r) " ©(OB + ^|e(?Xer) r"'e(M^,)) = f^ 
(aC(e,e(?/,^))Ve)-^^ = 5.' 
and the conditions on the contact boimdary 

r(z^,)e9/[o,i](^e^) J 
[De(Dtu^V) + Ce(u,^)-e(79^,)B + r|e(?/^,)p"'e(u^;)]z/± 

Kzllul] + (/!,)'( KM ) + [BeiBtuD + Ce(«^J,)-e(C)B 

+ r|e(.,^)r^e(.,^)).± + r'^(l+|K;] 1')^"'^] = 

i (D^ (e , e(^,^)) Ve l?c + ^(^er , e(^,^)) Ve Ire ) '^^ 

+ 77(K-i],z,^) [8(^)1 =0 
[3C(,?,^,e(«,^))VC];.± = -aiD^z.^ 
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To this aim, first of all we test equation (7.9a) by u^^. Thus, with elementary calculations, we find 

^ W k \\2 ^ i|A;||2 '^ll^'in 

+ .[ zJ,|Kl|^d5+/ (/^)'(K]).[uyd5 + r^/ {l+\lu',^]\Y'''\Hr]\'dS 

< C (^g\\u'^ ^llw"i.2(i2;Rd) + £'II"t ^llL2(0;Rrf) + ll-^THH'i-2(a:R'i)* + II /r II ffi/a (rN;R<i). ) 

+ 4r|Cn|e(u^)||i.(f,), (7.14) 

where d :— inf^gRd |^|=iIDC-C > 0, cf. (2.6), and where we have also used Kern's inequality in the form 

3CK^CKm yv eW^-^{n;R^) : \\v\\w^.^n:R^) < CK\\e{v)\\L2(n,R^..y (7.15) 

Also taking into account the monotonicity of the operator we have that the fourth, the fifth and the 

sixth term on the left-hand side of (7.14) are non- negative. Secondly, we test (7.11a) by z^^. Supposing 
that z^^ 7^ 2^^^, so that dC,i{{z'^^ — z^^^)/t) = — oi, we obtain with trivial calculations 

f / z^^\ [ulj I' d5 + r"||z,^J|i.(r,) + / riz^r)z^,dS < i2a^+ao)\\z^ALHTc) ■ (7.16) 

Note that the third term on the left-hand side of (7.16) is non- negative by monotonicity of the operator 
9/[o^i]. Finally, we test (7.9b) by -d^^, thus obtaining 

^lldli=(o, + K^jVePda; + 1^ r;(K\-i],z,^)[e(e)][e]d5 



< ^llC"'lli^(n) +Ii+l2+h+l4 (7.17) 



where we have used that 3C is positive definite (cf. with (5. Id)). As for the remaining terms Ij, j 
1, . . . , 4, we have 



- ^II^£tIIl2(o) + C* (||e(uJV)llL4(J2;E<'x<') + W'^i'^er ) II L4(0;R'' x''; 

1 fc 2 k k 1 

- ^ll^eTllL2(o) + ^^||tteT-||^yi.-,(j-j.jfd) -t- C||Mg.^ llH'i.T(0;Rd)7 (7-18) 

where we have used Holder inequality and the fact that 7 > 4. Furthermore, relying on (5.2), and setting 
Pi^ = 2uj / {lo — 1), we find 



r 



/2 = - / e(e) B: (e(«^;) - e{u^;')) d'^, dx 
" /o\rc 



<T^lieilW2)+C^ / \e{u^;r)-e{u';-' f\{^'i.r/''+l\dx 



16t' 
1 

87' 



< ^ll^?L||i2(f,) + C (||e(u^;)||^';^(,,^R.x.) + l|e(«.^, ')llip.(0;R^><.) + 1 



- ■^W'^srWhin) + ^^ll^eTll^vi.T(a;Rd) + C* (llitjr ^Ww^-^in-M") + 1) ' (7-19) 

where we have successively used Holder's and Young inequalities, and that 7 > p^j due to our assumption 
that 7 > max{4, -^jzj}- Besides, using that < z^^ < 1 and < z^'~^ < 1 a.e. in Fc, and the continuous 
embedding (4.1), we find 

^3 = -ai / B.z^^^^zllc^^izkc^s < C\\B,z^ALHra)\\^U\LHra) < pMr\\w^-^n)+ 
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where we choose pi > [Cp-^ being some constant depending on pi > 0) in such a way as to absorb 

^eTllvFi.2(n) 



I'^^erll wi.sfni t'^*^ left-hand side of (7.17). Finally, we have 



(7.20) 



in which we again choose a suitably small p2- Collecting (7.14)-(7.20), we readily conclude an estimate 
for |!u^^||H'i.T(n;K<j), llz^^llLoo(rc) and ||i?^Vll wL^cn)- 

Since € W^^'^{n\Ta), we have that e W^'^{n\Tc) is a legal test function for (7.9b). 

Hence, we use recursively that 'd^~^ > a.e. in ft (starting from the initial condition — "Oq > a.e. 
in n, cf. (5.4d)), the fact that G'; > a.e. in n and > a.e. in dn (cf. with (4.6c)-(4.6d)), and that 
— aiDtz^^ > a.e. in Fc, and finally the property [z?^'^]~6(i?gV) = a.e. in fi, due to the fact that 6 in 
non-decreasing (cf. (4.4)). Thus, we conclude that [^^^t-]^ = a.e. in f2, whence i?^^ > a.e. in 51. 

Finally, we prove (7.13) by adapting to the time-discrete setting a comparison argument from [19, 
Sect. 4.2.1]. Exploiting the fact that G^^ > a.e. in Jl, we deduce from (7.9b) that 

- div(aC(e,e(iiJ;))VO > ^|e(D,H,^)|' - C|e(^J,)p > -C'l^J.p in n\T^ (7.21) 

for any fc = 1, . . . , and some C" > independent of t and e, where d > is the positive-definiteness 
constant of D, cf. (2.6), and where the last inequality ensues from (4.4) and (5.1b). We compare (7.21) 
with the finite difference equation 

DtXfe = -C'lxfeP Vfc = l,...,i^„ (7.22) 

with C being the same constant as in (7.21). In fact, this is an implicit discretization of the ordinary- 
differential equations X + C"|xp = which, for x(0) = ho{9*) > with 9* from (5.8), gives a sub-solution 
of the (continuous) heat equation. This initial- value problem has the solution xit) = l/{C't + l/ho{9*)) 
so that, in particular, x{ ) > 1/{C'T + l/ho{9*)) > on [0,T]. Now we solve (7.22) recursively starting 
from the initial datum xo — ^o(^*) > 0. In this way we obtain an approximate solution to the mentioned 
initial- value problem which, for r — > 0, converges uniformly on the considered finite interval [0,T]. In 
particular, for r > sufficiently small, we may take for granted that, say, 

^-^^*-- C-T+l/U) + l ^" Vfc = l,...,AV. (7.23) 

For every fc = \,...,Kt we subtract (7.22) from (7.21) (the latter supplemented with the boundary 
conditions (7.10c), (7.11d)-(7.11e)), and we test the resulting inequality by — (i^g^. — Xk)~ ■ Thus, for all 
fc = 1 , . . . , Kt , 

iD,((e-xfe)^)' < -(e-x^)"D*(e-xfc) < c'dcp - ixfen(C - x^)" < o (7.24) 

in r2\rc; the latter inequality also due to the previously proved positivity ?9^^ > a.e. in fi. Sum- 
ming (7.24) over fc = 1, . . . , K-r^ we easily conclude that (i?^^(x) — Xfc)~ = for almost all x G il and for 
every fc = 0, . . . , Kt-, whence i?^^ > Xfc ^ X* > a.e. in fJ. This concludes the proof of (7.13). □ 

Approximate solutions. In accordance with notation (7.2), for all t > we shall denote by 

• Uer, v^^, "der, I^r' ^^id Zjt, the picccwisc coustant interpolants of the elements {it^^j^j^, {'9^r}k=i-' 
and {z^jfji; 

• by UjT-, "deT, and z^r, the related piecewise linear interpolants. 



We shall now state the weak formulations of (the boundary value problems for) equations (7.9a), 
(7.9b), in terms of the interpolants so far introduced by using "discrete test functions" . Indeed, one verifies 
that for every X^-tuples {w^jfji C lFp;^(n\Fc; M'') and {w^jfji C W^^'^{Vt\Vc), the approximate 
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solutions (u^r T i^^sT , Zi^r , u^T , "i^er T Zer) fulfil thc following: 
the discrete (weak) momentum balance equation 



e(4^) + Ce(uer) ~ Be(i?er) + T|e(ueT) r"^e(T4,^)) :e(TX,) dxdt 

{nz^rlnerj + (/iO' ( M ) + (l+ | M 1')"^'^'^) -{vr] dxdt 



Sc 

/ gitgri- — t)-Vt- dxdt + g / 'Ugr{T)-Vr{T) dx 
er Jn Jn 



~P iio^r-VriT) dx + / Fr'Vr dxdt + / f^-VrdSdt (7-25) 
Jn Jq JSn 

(where we have used the notation wo.r = "'^ ^"'^ = uq), which can be obtained from (7.9a), (7.10a), 
(7.10b), (7.11b), and (7.11c) by using a suitable discrete "by-part" summation formula (cf. [44, Formula 
(4.49)]); 

the discrete (weak) enthalpy equation 

/ i}^riT)wr{T)dx+ / aC(e(t4,^),^er)V^er-Vuvda;dt 
Jn Jq 

f]{ {u^-^ ,Z,r) [e(4r )] [Wr] dSdt ~ ^ ^^Wr dxdt 

^ ^De(iier):e(4r) + 6(4r)B:e(4r)^ da;dt + ^ GrWr dxdt 
2 

again obtained from (7.9b), (7.10c), (7.11c), (7. lid), and (7. lie) by the use of the summation (cf. [44, 
Formula (4.51)]); 

the discrete flow rule of the delamination parameter (cf. (7.11a)) 



Q 



2 

aikr ^'"""' dSdt+ / §oWr{T)dx+ I g^Wr dSdt , (7.26) 



9/(-oo,0] (4t(0) + f I I^feT(t)] 1^ + T"Ze^(t) + r{Zer{t)) - flQ - fli 9 
r{z,r{t)) e dlloM'^erit)) 



on Ec • (7.27) 



A-priori estimates. Like in [44, Lemma 4.1], we derive some further energetic information on the 
approximate solutions (see Lemma 7.6 later on) by recurring to an auxiliary minimization problem. 
With this aim, we first proceed to the validation of a suitable (strict) semiconvexity property of the 
stored energy. 

We further introduce the short-hand notation for the regularized stored energy 



^er{u,z):=. f flce{uy.eiu) + -\e{u)p)dx + II,{lul) 

+ 1 {^z\M\' + ^\z\' + jMluj\Y' + IloM^)-'^oz)dS. 



Lemma 7.5. Under the assumptions of Theorem 6.1, suppose further that the exponents fi G (4,5), 
a, /3 € (0,1) in (7.12) comply with 

a(M-2) + 2/?< (7.29) 



2 



Then, for every k > there exists > such that for all < t < the function on VFp^ (fi\rc; M'*) x 
L°°{Tc) given by 

/ \ f Be(u):e(w) 

(u, z) 1-^ $£^(7i, z) + / -= dx is strictly convex. (7.30) 

Jn 
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Proof. Following the calculations in [31], wo prove (7.30) by investigating the monotonicity of the multi- 
valued mapping 

To this goal, we have to estimate from below 

{d^eriui, Zi) - d^eTiu2, Z2), (ui - U2, Zi - Z2)) 

~ I (Ce(ui— M2) + T|e(ui)p^^e(ui) — T|e(u2)|'''~^e(u2)):e(ui— U2) da; 
Jn\Tc 

+ I [u2],zi,z2)d5. (7.31) 

Jtc 

As for the latter term, using the short-hand notation Si = fuil for i = 1,2, and r{z.i) E dI[Q_i]{zi) as 
in (7.11a), we can estimate the last term in (7.31) as 

L{si,S2,Zi,Z2) = K{ziSi-Z2S2)-isi-S2) + ( (/^ )' («! ) - (Z^)' (^2)) -(51-52) + ^ (zi -Z2) ( | Si p- | S2 P) 

+ (r(zi)-r(z2))(2i-^2) + r"|;^i-Z2p + r^" {il+\s,\^r^'-\si-il+\s2\'r^'-' S2)-is,-S2) (7.32) 

> KZi|si-S2p + |(21-Z2)(S1+3S2)-(S1-S2) + t"(zi-Z2)^ 
+ T^{{l + \s,\'r/'-'s,-{l + \s2\'r/^''s2)is,~S2) 

> Yi^i~z2r - ^\-si-s2ns,+3s2\' + r^{iiHsi\'r^'^'s,^a+\-s2\'r^'-'s2)isi-s2) 

t" 2 1 2 

>-^{zi-Z2) - S'K,r|si-S2| (7.33) 

for some positive constant — S'^.t- Indeed, the first inequality follows from the positivity (by monotonicity) 
of the second and fourth term on the right-hand side of (7.32), and from simple algebraic manipulations. 
So does the second inequality. To conclude the final inequality (7.33) for some constant Sk,,t > 
depending on k and r, we have used that (cf. [36, Lemma 5.2]) 

3C^>0 : (l + |si|2)^/2-l.Si-(l + |s2|2)'^/2"ls2)-(si-S2) > C^dsil^-' + |S2|^~')|S1-S2|' 

for all si, S2 G K'' and that (since > 4) 



k2 



yK,T>0 3 5,., >0 Vsi,S2eM'^ : c,,r^(|si [''-^ + |52|^-^)- — |si+3s2|' > -5„., • (7.34) 

OT 

Combining (7.31) with (7.33) the boundedness of the jump operator u i— > |u 1 from Wr^f (f^\rc;M'*) onto 
i^(rc;M''), as well as using Korn's inequality, we conclude that 

r" 2 

(d^er{ui, Zi) - 5$er(w2, 2:2), ("l " "2, ^1 " Z2)) > I kl "^2 1 1 ^2 (r^j -C'S'k.t || e(wi )-e(u2) || ^2 (fj-Rdx d) 

with the constant C depending on the positive-definiteness constant of C (cf. (2.6)), on the norm of 
the trace operator from W^^'^(r2\rc; M'') to L^(rc;M''), and on the constant in Korn's inequality (7.15). 
Finally, the key observation is that, for k > fixed, the constant S^.t in (7.34) has the following qualitative 
behaviour ^ 

Sk,t ^ — n-2 I 2/3 as T — > 0. 

Thus, using condition (7.29), it can be verified that for all k > there exists > such that for 
< r < T„ there holds CS^^r < d/-y/r; again d > is the positive-definiteness constant of B. This 
yields (7.30). ' □ 
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Lemma 7.6 (First a priori information). Under the assumptions of Theorem 6.1, for all g > and foi 

every e > there is such that for all < t < the approximate solutions {u!,T-,i!feT,ZeT,UeT,'&eT, z^t) 
fulfil the following "discrete mechanical energy" inequality 



^ ^kin(^0,T) + $eT('"o,r,2o) 



e(4r(s))B:e(uer(s))dx+ / Fr{s)-Uer{s)dx+ / f^{s)-Uer{s)dS)ds, 

as well as the following "discrete total energy'' inequality 



(7.35) 



+ 



( f Fr{syu,r{s)dx+ f (s) • 4r (s) d5 + [ Gr{s) dx + [ g,{s) ds] ds (7.36) 
\Jn JFn Jan J 

and also the "discrete semistability" for a. a. t e (0,r) (where is from (7.28)j 

^„{u,r(t),z,r{t)) <^,r{v^r{t),~z)+^{~z-zer{t)) for all S £ (T a) . (7.37) 

Proof. Let us now fix a solution (w^^, z^^, i?^'^) of Problem 7.1. Recall that such a triple exists thanks 
to Lemma 7.4. Let us consider an auxiliary minimization problem, namely 



minimize / gD^u^^-u + {l-y/T)Be(ptu^.^):e{u) 

- ' "^):e( ) + 9(z9^;)B:e(^)dx 



^/ Ci( ^ ^dS + <Per{u,z)- f F^-udx- f f^-udS 



subject to {u,z) £ W^2i^\^c''^'^) ^ L^i^c) ■ 



(7.38) 



By convexity of and coercivity (cf. the calculations developed in the proof of Lemma 7.4), it is 
immediate to check that the minimization problem (7.38) has a solution which we denote by {u^^,z^^); 
of course, it depends on the pair (u^^,i?^^) in general. Writing optimality conditions for (m^^,z^^) gives, 
for ah V e Vyji;7(17\rc;M'*), that 



) +Ce(ii^V) +T|e('S^^)|'' ^e(uj^)) :e(w) dx 



{u^j + mm^] ) + (l + I K^.] K?,] ) • H dS 



= / (l-V^)l])e(Dtu^;):e(w)-e(i5^JB:e(z;)+F^^-wdx+ / f!^-vdS 
Jn 

and for aU S £ L°°{rc) 



(7.39a) 



Ci(5)+Vz,^,-V 



dS" 



rc 



r"5^, + r(5^,) + - 



k 1 



dS* > 



^rc 



<■( 



dS*. (7.39b) 



Now, we test the difference of (7.9a) and (7.39a) by u^^ — u^^ and the difference of (7.11a) and (7.39b) by 
z^^ — z^^ and sum up the resulting relations. Using that the underlying potential, namely the functional 



(u,z)k>$,^(u,z)+ / Ciiz-z^-^)dS + 



De{u):e{u) 



dx. 
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is strictly convex on Wp;f(rj\rc;M'') X i°°(rc) by Lemma 7.5, we conclude that 



Then, the functional in (7.38) must have a bigger or equal value on {u 
which gives a discrete analog of (4.17), namely 



k-l _/£-l\ 



than on (u^^, z^^.) 



^,r{4r.-4r)+r I Ci{DtZ^r) dS + ^ f (2- )lD)e (D* 7.,^) :e (D^O dx 



£r 7 er 



r / e(C)B:e(Dtu^,) + F^-Btu^^.dx + t [ /^D^u.^dS (7.40) 



when also employing the algebraic inequality 

Upon summation over k, we conclude (7.35). 

Now, to get (7.36), we add to (7.40) the relation obtained testing (the weak formulation of) the 
boundary- value problem (7.9b, 7.10c, 7. lid, 7. lie) by r. Developing all calculations, one sees that, 
thanks to our carefully designed discretization, the fourth term on the left-hand side of (7.40) and the 
first dissipative/adiabatic term on the right-hand side of (7.9b) mutually cancel out. So do the third term 
on the right-hand side of (7.40) and the second right-hand-side term in (7.9b). Again upon summation 
over k, we arrive at (7.36). 

Finally, to check (7.37), it just suffices to realize that the functional minimized in (7.38) has a lower 
value in (u^^,z^^) than in (uj^,z) for any z € L°°(rc), which gives 



,fc-i 



dS. 



Then, by using that Ci is homogeneous degree 1 and thus satisfies the triangle inequality C,i{z—z^^ ) < 
Ci{S-z^^) + Ciiz^r-^T^)^ we find 

{ul, zl) < $er ~Z) + [ Cl {S^Z^r') ^ Cl {z^r~4r') dS < (l^, + / Cl {S-Z^r) dS . 

Being k ~ 1, . . . , Kr arbitrary, we conclude (7.37). □ 

Lemma 7.7 (A priori estimates). Under the assumptions of Theorem 6.1, there exist constants So and, 
for every 1 < r < such that for all g > 0, e > and for all Q < t < (V„ being as in 

Lemma 7.5), for all approximate solutions {JI^t , -d^r ,'ZeT , u^t , "^st t Zst) the following estimates hold 



l^'^IL~(0,T;Wj;'^(f2: 



||"eT||^i.2(o,T;H'p^''(a;K'')) 
o^'Hu. II 



<5o. 



We' 



^llL~(0,T;H'p'"'(n; 



< 



So 



■^llBV([0,T];Li(rc)) 



< Sr for anyl<r 



\L^{0,T;W^'^{n)) 



Li(0,T;lVi.'-'(f2)*) 



BV([0,T];H'p'"'(n;R<*)*) 



<5n 



(7.41a) 
(7.41b) 
(7.41c) 
(7.41d) 

(7.41e) 
(7.41f) 
(7.41g) 
(7.41h) 

(7.41i) 
(7.41j) 



where Sq and Sr neither depend on e nor on t. Estimates (7.41e), (7.41f), (7.41g), (7.41h) respectively 

hold for Zer, 'ZsT, l^er and der, CLS Well. 
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Proof. Some of the calculations wc shall develop hereafter arc analogous to the ones in the proof of [44, 
Prop. 4.2], to which we shall systematically refer. 

First of all, we use the "discrete total energy" balance (7.36). Indeed, on the one hand, by def- 
inition (7.28) of ^sT, the second term on the left-hand side of (7.36) is non-negative and, thanks to 
positive-definiteness of C and Korn's inequality (7.15), it provides a bound for ||t^T(i)||^i.2(Q.a<i) ^'^d for 
r|ji%i-(i)||^i,T(Q.][jd) uniformly w.r.t. t g [0,T]. Further, being > a.e. in Q thanks to (7.13), the 

third term on the left-hand side of (7.36) estimates ||i?er |!l°°(o.T;L1(0))- To estimate the right-hand side 
of (7.36), we employ the discrete "by-part" summation [44, Formula (4.51)], to the effect that 

iAt) r _ r _ r _ r^At) r . 

/,(s)-4r(s)d5ds- / U{tyu,rit)dS^ (t) -Uo.r d^- / / (s) -U., (s) d^ds 



< 



P3||'"er(t)|j^i.2(s2;Rd) + C'pa (ll"0,T||^i,2([2.Rd) + II /r II L~(o,T;Li/3(rN))) 



C / ||/r(s)||LV3(rN)ll^r(s)||M/i,2(t,.Rd)ds, (7.42) 







where inequality (7.42) is also due to the continuous embedding (4.1) and ps is chosen in such a way 
as to absorb the first term on the right-hand side into the term |ji'eT(t)||^i,2(f2.Rd-) on the left-hand side 
of (7.36). Furthermore, in the case g > we estimate the fourth term on the right-hand side of (7.36) by 

Fr{s)-iMr{s)dxds < / ||Fr(s)||L2(0;Rd) ||wer(s)||L2(n;R'') ds . (7.43) 
Jo 

We then combine (7.36), (7.43), and (7.42), and use (7.5)-(7.6b) for TV, G^, and g^. Applying 

the Gronwall Lemma, we conclude estimates (7.41a), (7.41c), (7.41d), (7.41e), and (7.41g) (the estimates 
for z^r and following from the bounds for z^t- and i^^t-, and from (7.3)). In the case g ~ 0, the only 
change in the above calculations is that, under the additional assumption (5.5a), we estimate the fourth 
term on the right-hand side of (7.36) by use of the aforementioned discrete by-part summation formula. 
Namely, on account of the Sobolev embedding (4.1) 

/ / Ft-{s)-u^t{s) dxds ^ / Fr{t)-u^rit) dx — / Fr{T)-uo^r dx — / / Fr{s)-ii^^{s) dxds 
Jo Jn Jn Jn ' Jt Jn 

< P4||'fer(i)||^l,2(f2.Rd) + (^ll"0,r||^l,2(f^.Rd) + 1 1 || ^oo (o,T;L6/5(n;M<' )) 

f-uit) . _ (7.44) 

+ C ||F^(s)||i6/5(Q.Rd)||Uer(s)||M/i,2(0;Rd)ds, 







where again the positive constant is such that the first term on the right-hand side of (7.44) is controlled 
by II ^T(0llvi/i,2(Q.Rd) on the left-hand side of (7.36). 

Secondly, again arguing as for [44, Prop. 4.2], we make use of the technique by Boccardo & Gal- 
LOUET [7], with the simplification devised in [18]. Hence, we test the heat equation (7.26) by H^-der), 
where 11 : [0, +oo) — > [0, 1] is the map 

w^U{w) ^1- <?>0; 
(l+w)« 

note that Tl^-der) G W^''^{rt\Tc), because 11 is Lipschitz continuous. With the same calculations as in [44], 
taking into account (5. Id) we find 

<;k / —I^Jll!— dxdt < / 3C(e(l2eOJer)V4r-Vn(4r)dxdt 

+ f Vilier,-^er) [e(^er)] {lli^er)} dSdt + f n(4r (T, •)) dx 
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< / n(^?o)da- + ||a|Ui(Q) + l|ff,llLi(s)) 
Jn 

+ C(||ee(4.):e(4r)|Ui(Q) + |ie(4.)IB:e(^}er)|lLi(Q) + |lCi(4r)|!Li(Sc)) (7.45) 

where 11 is the primitive function of 11 such that 11(0) = 0. Note that inequahty (7.45) follows from the 
fact that ?7(Mg^,^gT-)[[8('i?eT-)l Iln(t9er)l > a.e. in Ec (by the positivity of rj and the monotonicity of Q and 
n), from the "discrete chain rule" [44, Formula (4.30)] for ft, and from < U^der) < 1- Combining (7.45) 
with the Gagliardo-Nirenberg inequality, we find, for all 1 < r < {d+2)/{d+l), that 

ll^'^-lll^(Q;R^) ^ ^r{l + ||©e(4r):e(4r)||Li(Q)+ m%r)M:e{{l,r)\\ LHQ) + 1 1 Cl (^er ) |1 (Sc))) (7-46) 

for some positive constant Cr, depending on r and also on the function 77, cf. (5.1e). 

Then, we multiply (7.46) by a constant > and add it to (7.35) (in which we set t = T). Now, 
by positive-definiteness of C, the third term on the left-hand side of (7.35) is bounded from below by 
d/2||e(?4r)||^2(Q.Rdxd) J whereas the fourth term controls ||Ci(iET)||Li(Ec) ■ Thus, we choose small enough 
in such a way to absorb the first and the third term on the right-hand side of (7.46) into the left-hand 
side of (7.35). Hence, we find 

^l|e(?4r)|li2(Q.Rdxd) + (l-p5)||Cl(2£r)||Li(Sc) + P5 1 1 V^r 1 1 ^.(g.^d, < T^^j^ (uo,r ) + $£r (wo.r , ^O.r ) 
T 



Fr-i^rdxdt+ / /,-4rCl5dt+(p5a + l)|le(l?,,)Be(Uer)|lLi(Q)- (7.47) 

Jq JTii 

The first two terms on the right-hand side of (7.47) are estimated in view of (7.5)-(7.6b) and (7.7), 
whereas, taking into account the Sobolev embedding (4.1) and Korn's inequality (7.15), we have 

2 , ° / M|2 



Fr-ikTdxAt<C I \\Fris)\\^e/^n.xd) + — ||e(^}er)|lL2(f2.Rdxd)dt, 







T 



d 



U-V^r dSdt <C ll/r(s)lli4/3(r„;Rd) + l|e('4r)|li2(0;K<ixd) dt] 

here again d > is the positive-definiteness constant of D. To estimate the last summand, we use the 
positive-definiteness of C and (5.2), finding 

(p5a+l)||e(^,0B:e(ii,0llLi(Q) < P6|B| ||e(?i,.)lli2(Q.Rdxd)) + C,J|e(4r)|!i2(Q) 

< P6|D| ||e(4r)||i2(Q^Rdxd)) + Cp, (||4r||l{/.(Q)+ 1) (7.48) 

in which we choose the positive constant pg small enough, again to absorb the first term on the right-hand 
side of (7.48) into the left-hand side of (7.47). In order to estimate IjiJer |!l2/"(q)j we again employ the 
Gagliardo-Nirenberg inequality. Indeed, with the same calculations as throughout [44, Formulae (4.39)- 
(4.43)], and relying on the restriction of w in (5.1b) and on the bound for ||'i?eT||L~(o,T;Li(f2))j '^e conclude 

^ mfl.:.^^)) < pr f PMLiu-M^) + (7.49) 







for a suitably small py > 0. Then, we plug (7.49) into (7.48), and the latter into (7.47), and choose 
Pq in such a way as to absorb II ^^^er H^r^g.jjd) into the left-hand side of (7.47). Thus, we conclude 

estimate (7.41b), as well as an estimate for ||C(^er)||Li(Sc) (yielding (7.41f)), and a bound for S/der in 
L^iQ; R'^). Combining the latter information with the estimate for in L°°{0, T; L'^{i})), we infer (7.41h) 
(the estimate for •d^r due to the bound for -der and to (7.3)). As a by-product of the above calculations, 
we find 

IIAerllLi(Q) < C, with A,r ^—T^V)e{^ry-e{^r) + &{^er)B:e{{,^r) ■ (7.50) 
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For later convenience, we also remark that (7.41h) yields 

(If 

l|e(4r)||L^-(o,T;L-.(rc)) <C for all 1 < 9 < j— , (7.51) 

where we have also used the continuous embedding W^'^{il.) C L'^{rc) for q ranging in the above- 
mentioned index interval, as well as the growth restriction (5.2) imposed on Q. 
To prove (7.41i), we argue by comparison in (7.26), to the effect that 

ll'^e^llLi(o,T;Wi.'-'(n)*) = sup {h + Ie+ h + h) , where 

= (0,T:H'l-'"' (SI)) — ^ 

h= ^tW < ||Aer||Li(Q)||w||L-(Q) < C\\w\\ r^) 

JQ 

thanks to (7.50) and the continuous embedding W^'"^ (O) C L°°{Vt) (since r' > d + 2), while 



3C(e(UeT-), ^?er)Vl?e^-Vw < CgcH V??er || L^IQ) || Vti;|| (qj < C\\ 

I 

due to (5.3) and (7.41g). Further, 

/r = - ^ ((zer [e(4.)] H + ai4.^^^^4^) '^^'^^ 

< (ll%r|U°=(Ec) l|Q('^er)||L-'-(0,T;L-<J(rc)) + II Cl (4t) || Li(Ec)) lk||L~(Ec) < L-^ {0,T;W^y {n)) 

thanks to (7.41e), (7.51), and the continuous embedding W^'^ (17) C L°°{Tc)^ and, finally, 

/g = / GrWAxdt+ / g^wdS'dt < (llGrllii(Q) + ||w||i=o(o,T;H'i>'-'(0)) ■ 

Collecting the above calculations, we conclude (7.41i). 

Finally, for (7.41j) we use that u^t is a measure on [0,T], supported at the jumps of ii^T^ and we 
estimate the norm g|| iisr |Im(o T-w-'^-"'(r!\rc-R'')*)' '^^ere M(0, T; M^p^'''(ri\rc; M'')*) denotes the space of 

Radon measures on [0, T] with values in iyp^"'^(ri\rc; M'*)*, arguing by comparison in (7.25); see the proof 
of [44, Prop. 4.2], where similar calculations were carried out. □ 



8. Limit passage with r — >■ and proof of Theorem 6.1 

Throughout this section, we shall keep e > fixed, and let t — > 0. We shall develop a proof of the 
passage to the limit unifying the cases g > Q and g — Q. 

Step 0: selection of convergent subsequences. First of all, it follows from estimates (7.41b), 
(7.41c), and (7.41j), from the Banach selection principle, and from the Aubin-Lions theorem (see, e.g., 
[47, Thm. 5, Cor. 4] and [41, Cor. 7.9] for the generalization to the case of time derivatives as measures), 
that there exist a (not relabeled) sequence t and a limit function G VF^'^(0, T; Wj';'^(f2\rc; K'')) 
such that the following weak, weak*, and strong convergences hold as t — > 0: 

uer^u, u-vW^-''{Q,T-Wl'^{n\Ta-X)), (8.1a) 
iter^ue inC°([0,T];Wp^-^''(ri\rc;K'')) Vee (0,1], (8.1b) 
ifg>0, Uer^ue inW^^°^{{),T;L^{Vl-m.'^)). (8.1c) 

In the case ^? > we also have u^r in W^^^{Q,T]WI-^''^{VL\Tc;M.'^)) n W^''i{Q,T]L^{n]W^)) for all 
e G (0, 1] and 1 < q < oo. Furthermore, estimate (7.41j) and a generalization of Helly's principle (see [4] 
as wen as [35, Thm. 6.1]) yield that iie £ BV([0, T]; Vrp;'^(rj\rc; M'')*) and, in addition, iier{t)-^ Ue{t) in 
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Wj!;^'''(r2\rc; M'')* for all t G [0,r]. By virtue of estimate (7.41c) and of a trivial compactness argument, 
this pointwise weak convergence improves to 

iierit)-^ iieit) in M'') for all t G [0, T], in the case g > 0. (8. Id) 

Combining (8.1a) and (8.1b) with the general inequality (7.4), we conclude that, up to the extraction of 
a further subsequence, for all e e (0, 1], 

%r^u, mL^{0,T;W^-^{n\ra;R'')), u,r ^ u, in L°-{0,T;W^-''^{n\ra;R'')), 

lkr{t) u,{t) in W^]^^''2(l}\rc;IR'')) for all t € [0,T], ' ^ 

the latter pointwise convergence due to (8.1b) and (7.4). 

With the aforementioned compactness results, we deduce from estimates (7.41e), and (7.41f) that 
there exists a function S L°°(Tjc) ^ BV([0,T]; Z), {Z being some reflexive space such that L^{Tc) C Z 
with a continuous embedding, for example Z = W^'^~^'^{Tc)* for some e > 0), such that (possibly along 
a subsequence) 

z,r^z, ini°°(Sc), (8.2a) 

and, again by [35, Thm. 6.1, Prop. 6.2], lerit)^ Ze{t) in Z for all t E [0,T]. In view of (7.41e), we indeed 
have pointwise weak* convergence in L°°{Tc), i.e. 

z^rit) ^ z^{t) in L°°(rc) for aU t e [0,r]. (8.2b) 

Finally, arguing as in the proof of [35, Thm. 6.1] with Helly's selection principle and taking into 
account (8.2b), we conclude that for all < s < t < T 

Var3j(z,;[s,i]) < lim / / Ci (4r (r)) d^dr , (8.2c) 

which ultimately yields that there exists S'q > such that 

eBV([0,T];Li(rc)) and ||z,||Bv([o,T];Li(rc)) < for aU e > 0. (8.2d) 

Thirdly, by the same tokens we conclude from estimates (7.41g), (7.41h), and (7.41i) that there exists 
i?e e L''(0,r;VFi''-(rj\rc))nBV([0,T];M^i''^'(rj\rc)*) such that 

ini'^(0,T;W^i''-(O\rc)), (8.3a) 
4r,4r^i?e m V {0,T;W^'''-'{n\rc)) D L" {0,T; 1^^)) , (8.3b) 

for all e E (0, r — 1] and I < q < oo, as well as 

^erit), i?er(i)^i^e(i) m W^'''' {n\Tc)* for all t e [0,T]. (8.3c) 

Notice that, under condition (5.8) on ^o, convergence (8.3b) and (7.13) yield (6.8). It also follows from [35, 
Thm. 6.11 that Varii. II , (tJ^; [s, tl) < lim^_j.n VariLii , (iJ^t-; [s, tl) , with VariLii , denoting 

the total variation w.r.t. the norm || • This entails that 

ll^ellBV([o,T];Wi.'-'(n)-) < -^o for all £ > 0. (8.3d) 
For later purposes, we also point out that, in view of estimate (7.41g) and of (8.3b), there holds 

¥e\\L'^{o,T:LHn)) < 5o for all e > 0, (8.3e) 
Sq being the same constant as in estimates (7.41). 
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Besides, (7.41d) yields that 

T\\\ei%r)r-^e{u,r)\\L-,/(-y-i)(^Q.u^.^^<SoT^^^ ^0 as T^O. (8.4a) 
In view of (4.1) and the second of (8.1e), it is not difficult to verify that, for all e G (0, 3], 



M inL-(0,r;L4-(rc;M'^)), 

lu,{t)} in L4-<^(rc;K'^) for all t G [0,r]. 



.4b) 



Furthermore, using that (I^)' is given by (6.3), and recalling (6.4), from (7.41b) we easily infer that 

35i = 5i(e)>0 Vr>0: )|L^(o^^^^.(p^^„.„ < ^i, (8.4c) 

with 5'i(e) oo as £ — )• 0; more specifically, due to (6.4) we have Si{e) ~ ^'{l/y/e). Combining (8.4b) 
with the strong- weak closedness of the graph of the operator (Z^)', up to the extraction of a further 
subsequence we find that 

(/i)'(M) (/iO'(M) in L^iO,T;L\Ta;R'')). (8.4d) 

Moreover, using that |(l + IHuIttII^) ^ ^II^%rl| < 27^2^|u:j^^| | _|_ "') a.e. in Sc, as well as esti- 

mate (7.41b), one sees that the sequence {(l-h| |?^rl P)^"ni^rrl } is bounded in L°°(0, T; i^/^^-i) (Fc; M'*)). 
Thus, 

T^/\l + \lTkrj\')^''lTkrj ^0 inL-(0,T;L4/(/^-i)(r,;M'^)). (8.4e) 

Next, let us point out that, in the case the space dimension is d = 3, (7.41h) holds for all 1 < r < 5/4, 
so that (8.3b) yields by interpolation 

4r ^ i?e in L^^^'^^'iQ) for aU e G (O, f ] . (8.5a) 

In particular, 0(i?j^) -> 8(i?e) a.e. in Q. Combining this information with (5.2) (note that, by (5.1b), 
(jj > I for d = 3), it is immediate to deduce from (8.5a) that, for example, 

&(d,r) e(i?e) in L^iQ). (8.5b) 
Furthermore, using standard trace theorems we also deduce from (8.3b) that for all e G (0, |] 

^,+ |rc^^?lrc and ^;^\r^^^7\ra in V{Q,T; L'''/'-%Ta)), 
so that, again by (5.2), for d ~ 3, using that cj > 6/5, we conclude that 

[6(40] ^ [0(^e)] in L'-(0,r;ii2/7-.(re)) Ve G (O, f ] . (8.5c) 

Similar calculations leading to (8.5b) and (8.5c) can be performed in the case d = 2. 
In the end, we are now going to show that 

$e(we(i),^e(i)) < lim inf $e^(i%^(t), (t)) for aU t G [0,T]. (8.6) 

r— >0 

Indeed, taking into account (8.1e) it is not difficult to deduce that for all t G [0,r] 

liminf / ^Ce{T^r{t)):e{%r{t))dx> j ^Ce{u^{t)):e{u^{t)) dx. 
Combining (8.2b) with (8.4b), we have 

liminf / '^z,rit)\luer{t)}f dS > f ^z,{t)\lu,{t)}\^dS for ah < G [0, T] . 
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7b) 



Besides, taking into account that If^ is lower semicontinuous on L^(rc;ffi'') (cf. (6.1)), wc immediately 
conclude 

lim inf If, ( lu,r (t )]) > ( [ue (t )] ) for all t G [0, T] . 
Collecting the above inequalities and also relying on (8.2b), we infer (8.6). 

Step 1: passage to the limit in the momentum equation. As a first step, we shall take the 
limit as r — > of the discrete momentum equation (7.25) and of the discrete heat equation (7.26) with 
more regular test functions, which, for technical reasons, we shall need to approximate carefully. More 
precisely, for the momentum balance equation (6.7) we shall use test functions 

V e L^{0,T;W^'^^"{n\rc]R'^)) n W^'\0,T;L'^{n;R'^)) for some v > 0, (8.7a) 

and we shall approximate them with discrete approximations {w^:}, such that the related piecewise con- 
stant and linear interpolants fulfil as t — > 0, 

Vr^v mW^'\Q,T;L'^{n;R'^)), 
Vr-^v in L'^{G,T-W^'^+^'{n\Tc.-R'^) for some u > 0, 

Ti/7e(t7^)^0 in LT(Q; «'*><'*), 

'^''^^l|Wr||Li(0,T;Li/(5-M)(rc)) - ^ ■ 

Now, combining (8.1a) and (8.5b) with the second of (8.7b), and (8.4a) with the third of (8.7b), we pass 
to the limit as r in the first integral term on the left-hand side of (7.25). Secondly. (8.2a) and (8.4b) 
yield 

KZer {ver} ^ ^z, {u,} in i°°(0, T; L^''{Ta)) for all e e [0, 3), 

which we combine with the second of (8.7b). Also taking into account (8.4d) and (8.4e), together with 
the fourth of (8.7b), we take the limit of the second integral term on the left-hand side of (7.25). In 
the case p > 0, we take the limit of the third and fourth terms on the left-hand side, and of the first 
term on the right-hand side of (7.25) by means of (8.1c) (combined with the first of (8.7b)), of (8. Id), 
and of (7.7). Finally, using (7.5a) and (7.5b) we handle the second and third right-hand-side terms 
in (7.25). We thus conclude that the triple [u^^z^^-d^) fulfils equation (6.7), first with test functions 
as in (8.7a) and ultimately, by a density argument, with test functions v G i^(0, T; W]i'^(51\rc; M'^)) H 

Step 2: passage to the limit in the semistability condition. We consider a subset 3Nf C (0, T) of full 
measure such that for alH e ?^ the approximate stability condition (7.37) holds independently of r 0. 
Then we fix t e 3\f and z E L°°(rc). We may suppose without loss of generality that 3?(z — ^^(t)) < +oo, 
hence 

z{x)<z^{t,x) fora.a. xSFc. (8.8) 
We then construct the following recovery sequence 

{_ z(x) 

Zer {t, x) where z^{t,x) >0, 
Ze{t,X) (8.9) 

where z^it, x) = 0. 

Now, using (8.8) and (8.2b) one immediately sees that 

Ser{-,t) <z,r{-,t) a.c. iuTc, Ser{t)^SmL°°{Tc). (8.10) 
Plugging z^r in (7.37), we find 

< lim sup Ueri%r{t), S^r{t))+^{Ser{t) -Z,r{t))-'^>er{'Uer{t),Z^rit))] 

= limSUp^ (y (l^er(t)|'-|2er(t)|') + ^| [wer(i)] |' (^er(t) " %r(t)) - (co+ai ) (^er (t)" %r (t)) ) dS 

= ^e{u,{t), ~z{t))+'R{~z{t) - z,(t))-$,("s(t), z,{t)) , (8.11) 
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where the second inequality ensues from (8.10) and (8.4b). 

Step 3: passage to the limit in the mechanical and total energy inequalities. Using (8.1a), 
(8. Id), (8.2c), and (8.6), we pass to the hmit in the left-hand side of the discrete mechanical energy 
inequality (7.35) by weak lower semicontinuity. To take the limit of the right-hand side, we employ (7.7), 
the weak convergence (8.1a) and the strong convergence (8.5b). which yield 

e(4r)B:e(?4r) e(z9e)B:e(u^) weakly in ^^(Q). (8.12) 

Also using (7.5a)-(7.5b), we conclude that the triple (u^, Ze,'!?e) complies for all t G [0,r] with 



Jn 



+ $e(uo,^o)+ f ( f e{Ms))»A'^e{s)) dx + f F{s)-ueis)dx+ f f {s)-Ue{s) ds) ds . (8, 
Jo \Jn Jn Jr^ J 



13) 



We also pass to the limit in the discrete total energy inequality (7.36). Indeed, one tackles the left- 
hand side by the above-mentioned lower-semicontinuity arguments (also using (8.3c)), and passes to the 
limit in the right-hand side by convergences (7.7) and (7.5)-(7.6b). Thus, the total energy inequality for 
the e-approximate problem holds for all t G [0, T]. 

Step 4: mechanical energy equality. Like in [44], we prove that, in the limit, the mechanical energy 
inequality (8.13) in fact holds as an equality, obtaining 

TL{'^eit))+^e{ueit),zS))+ f f ne{u,):e{u,) dxdt + Yi,roiiz,;[0,t]) = T^.Juo) + ^^Uo, Zo) 

(8.14) 



Jn 



[ ( [ Q{i)e)V>:e{ue) dx + [ F-Uedx+ [ f{s)-Ue{s)ds]ds foralHe[0,T] 
Jo \Jn Jn Jfn / 



To this aim, we develop the same calculations as throughout [44, Formulae (4.69)-(4.76)]. The first step 
of the argument is a sophisticated trick based on the previously proved semistability condition (see also, 
e.g., [15, 24] for the use of such a technique in a rate-independent context), which allows us to prove the 
following inequality for all t € [0, T] 



$,(M,(i),z,(t)) 

t t (8-15) 

= // Ce{u,):e{u,)dxds+ f [ {Kz,[u,j -{u,} + (li.)' {lu,jy[u,j) dSds , 
Jo Jn\rc Jo JTc 

where {^e)'u denotes the partial Gateaux-derivative with respect to u of the functional : W^''^{Q\Tc) x 
L°°{Tc) and the equality follows from the definition (6.6a) of $£. The second step consists of testing 

of the momentum balance equation (6.7) by u^. In the case g = 0, e L^(0, T; VFp^^(f2\rc; M'')) is an 
admissible test function for (6.7). In the case > 0, the test by may be performed after proving that 
u, e L^{Q,T-Wl''^{n\Tc\W^)*), cf. Remark 5.3. In fact, a comparison argument in (6.7) readily yields 
that it£ e L^(0,T; W-^'_^{fl\Tc',R'^)*)- Choosing ii^ as a test function in (6.7) and integrating on (0,t) for 
all t e [0,T] leads, after an integration by parts, to 

Q 
2 



u^{t)\ dx + / I}e(u£):e{ue) dxds + / / Ce{ui.):e(ue) dxds 
n Jo Jn Jo Jn 



[KZ. 

JFc 

t 



W +(/!,)'( KD-W) dsds 

= 1 / |uo|^da;+ f ( f e{de)M:e{ue) dx + f F-Uedx+ f f{s)-Ue{s) ds] ds . (8.16) 
2 Jo Jo \Jn Jn JFn / 

Combining (8.15) with (8.16), we obtain the reverse inequality in (8.13) and thus conclude (8.14). 
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Step 5: passage to the limit in the enthalpy equation. First of all, wc observe that the following 
chain of inequalities holds for all t e [0, T]: 



Var3j(ze; [0,t]) + / / De(Me):e(ue) dxdt 
Jo Jn 



<liniinf / / Ci (4r) dS'di + / / De(iier):e(4r) dsdt 

< liniSUpTifij^(uo,r) + $£r(uo,r,2o) - T^ini'^rit)) " (l%T (*) , ^er (*) ) 

/ e(4r)B:e(?4r) + / F,-4rdx+ / 7,-4rd5) dt 
Jn Jn JTn / 



< T^kin("o) + $eK, 2o) - T^iniMt)) - ^e{ue{t), z,(t)) + / ( / e{d,)M:e{u,) + F-u, dx 



n 



+ / /-Ued^) dt = VarK(ze; [0,t]) + /" j De(ue):e(ue) dxdi. (8.17) 
Jfn / Jo Jn 

Indeed, the first inequality ensues from (8.1a) and (8.2c), the second one from the discrete mechanical 
energy inequality (7.35), the third one from (7.7), (8. Id), (8.6), (8.12), and from (7.5a)-(7.5b), cf. also 
Step 3. Finally, the last equality ensues from (8.14). Thus, all of the above inequalities turn out to hold 
with an equality sign. By a standard liminf /limsup argument, this entails that 

lim.^oT^J5.(t))_= TLM)), \ fo, , ^ [o,T], (8.18) 

lmv^O*er(Uer(t),Zer(0) = J 

as well as 

De(4^):e(4^) Be{u^):e{ue) strongly in L\Q). (8.19) 
Furthermore, arguing as in [44], from (8.17) holding as an equality we conclude that 

Ci(4r) S.f"^ in measure on E^T, (8.20) 

with ^1"''^ being the measure introduced in (4.15). 

We are now in the position of taking the limit of (7.26), where we shall use test functions 

w eC°{0,T;W^^'''+'{n\Tc))nW^^'''{Q,T;L'''{n)) for some <r > 0, (8.21) 

and we shall approximate them with discrete approximations {w^}, such that, r — > 0, the related inter- 
polants fulfil as ttv ^ w in C°(0, T; W^-'"'+' {n\Tc)) for some <r > 0, and ^ w in W^'''' {0, T; L''\n)). 
Then, we pass to the limit in the first integral term on the left-hand side by exploiting (8.3c) and 
the aforementioned convergence for the test functions Wer- To deal with the second term we observe 
that, due to (8Je), to (8.3b), and to the boundedness of the function % : M^^'' x M ^ M^^'', there 
holds 3C(e(i%7-), ??jT-) — ^ 'X{e{u^), d^) in L'^{Q) for all 1 < g < oo, which we combine with the weak conver- 
gence (8.3a) for and with the convergence for u%r. It follows from (5.1e) and from convergences (8.2a) 
and (8.4b) that rjHu^^j ,Zer) ^ r]{lu,j,z^) in L°°{0,T; L^~'{Tc)) for all e e (0,2], which we exploit with 
(8.5c) to take the limit of the third integral term. The passage to the limit in the fourth term results from 
(8.3a) and the convergence for i)gr- As for the right-hand side of (7.26), to deal with the first integral 
term we exploit (8.19) and the convergence for Wer, which in particular yields Wer — > w in C*'(Q). Relying 
on this convergence and on (8.20), we also infer 

hm (- /^^«r4. ^-'^-f aSdt) . £^^^^^4^ C"(d^dt). 

Finally, employing (7.6b), one takes the limit of the last three terms on the right-hand side of (7.26), 
thus finding that the triple (u^, Zg, 'd^) fulfils the weak formulation (4.14) of the enthalpy equation for all 
test functions as in (8.21). Again by a density argument, we conclude that (ue,Ze,?9e) fulfil (4.14) with 
test functions w £ C°{[0,T];W^^''' {n\Tc)) nW^-^'' {0,T; L''' (n)). 

This concludes the proof that solves (Me,^^,??^) the approximate problem, i.e. Theorem 6.1. □ 
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Remark 8.1 (Strong convergence). Let us observe that, in the case p > 0, by a classical argument 
based on a Korn-type inequality and the uniform convexity of the space W^''^{0,T;W^'^{^l\rc',^'^)), 
convergence (8.19) joint with (8. Id) allows us to conclude that u^r -> u, in W^'^{0,T; W^'^{fl\Tc]^'^)). 
Likewise, ii g > (8.18) gives {t^r{t) Ue{t) in L'^{n;R'^) for ah t e [0,T]. 

Remark 8.2 (Numerics). We point out that our method of proof may yield some strategy for numerical 
analysis after making a spatial discretization, although the non- variational structure of (7.9) (preserved 
if discretized in space) would still require some iterative procedure for numerical solution, cf. [5]. 

9. Limit passage with £ — and proof of Theorem 5.1 

In passing to the limit in the e-approximate problem as e ^ 0, we shall follow the steps of the proof 
of Theorem 6.1. Thus, we shall sketch most of the arguments, referring to the detailed calculations 
developed in Section 8, and dwell with some detail only on the passages to the limit as e — >■ in the 
momentum equation, and on the proof of the mechanical energy equality. 

Step 0: a-priori estimates and compactness. The sequence (u^, z^, i?e)e inherits the a priori esti- 
mates of Lemma 7.7, i.e. now 



ll'"'^llvyl-2(0,T;Wr^^(O;R<i)) + ll"^llwl.~(0,T;L2(O;R<i)) + £'11 "ellL2(0,T;lV^^^(n;R<i)*) (9.1a) 

ll^=IL°°(Sc) ll^'^llBV([0,T];Li(rc)) - ^ ' (9.1b) 

ll'^'^IL°°(o,T;Li(f2)) ~^ II'^'^IIl'-(o,T;W1.'-(o)) ll'^'^ IIbviio.tIiWI---' (o)*) - any 1 < r < , (9.1c) 

sup ^eiUe{t),Z,{t)) < S. (9. Id) 

te[o,T] 



for some 5 > and S'^ > depending on 1 < r < ^ij. Indeed, the first two estimates in (9.1a), the 
first of (9.1b), and the second of (9.1c) respectively follow from (7.41b), (7.41c), (7.41e), and (7.41h) 
via lower-semicontinuity. The second of (9.1b), and (9.1c) have been proved throughout Section 8, 
cf. with (8. 2d), (8.3d), and (8.3e). The third of estimates (9.1a) follows by testing (6.7) by functions 
V e L^{0,T-W^'^{n\Tc;M'^))riW'^'^{0,T]L^{n;R'^)) and taking into account (9.1a), the second of (9.1c). 
Finally, (9. Id) is a direct consequence of the total energy inequality. 

By the Banach and the Helly selection principles, there is a subsequence (for simplicity, denoted by 
the same indexes) and (u, z, t?) such that the following convergences hold: 



Ue^U 




(9.2a) 


Uf, ^ u 


xnC\%T]-Wl-^-\n\T^-W')) Vee (0,1], 


(9.2b) 


if £»> 0, 


u,^u in W^i'°°(0, T; L^{n\Tc; M'')) n ^^^'^(O, T; Wl'^{n\Ta] M'')*), 


(9.2c) 


if ^? > 0, 


u,^u in W'^'^{Q,T]W^'^'\a\Tc;«^)) n iyi'«(0,r;L2(f];K'i)) VO<e<l, l<g<oo. 


(9. 2d) 


if £»> 0, 


Ue{t)^u{t) in L^{n;M.'^) for aU t G [0,r], 


(9.2c) 


Ze ^ Z 


in i°"(Sc), 


(9.2f) 


Ze{t) ^ ■ 


z{t) in L°°(rc) for ah t e [0,T], 


(9.2g) 




in L'-(0,T;Wi--(f]\rc)), 


(9.2h) 




in L''(0,T;W^i^'-^(l]\rc))nL«(0,T; 2.1(17)) VeG (0,r-l], \<q<oo, 


(9.2i) 


§e{t)^^{t) in W^'''' {n\Tc)* for all t e [0,T], 


(9.2j) 


hem - 


lu{t)\ in L^-^V^-W') for any t € [0,7], J 


(9.2k) 




in L\Q), 


(9.21) 




^ [e(^)] inL'-(0,T;L3/2(r^)). 


(9.2m) 
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Convergences (9.2) can be deduced from estimates (9.1) arguing in the very same way as through- 
out (8.1a)-(8.5c) in Section 8. 

Step 1: passage to the limit in the momentum equation. First of ah, notice that (9. Id), (9.2k), 
and (6.2) yield S* > liminf^^^o /|<(lwe(t)]) > /^([[^(i)!) for aU t € [0,r], whence (4.10). We now exploit 
(9.2a-f,k-l) to pass to the limit in (6.7) with g>0, tested by ?; - u^, for any v in L^^q, T; W^'^{n\rc; M'')) 
(and in W^-^{0,T; L'^{n;R'^)) if g > 0), with [w] h on Sc, i-e. for any test function in (4.11). Using 
the mentioned convergences, we find 

limsup / (40'( dSdt = limsup / guo-{u,iO)-viO)) ~ gu,{Ty{u,iT)-v{T)) dx 

+ / fiu,-v)dSdt- f Kz,lu,j-lu,-vj dSdt 

+ / F-(u^-v)- (©e(weJ+Ce(we)-Be(i?eJ):e(we-w) + eiwV(we-v)da;dt 
JQ 

< I guo-{uo-v{0)) - gu{T)iu{T)-v{T))dx + [ fiu-v)dSdt- [ kz [u] • d5dt 
Jn JSn J'Sc 

+ / F-{u-v) - {l]ie{u)+Ce{u)-Me{-d)):e{u~v) + gu-{u-v)dxdt. (9.3) 

On the other hand, recalling formula (6.3) for the Yosida regularization (Z^)', we see that 

"/(40'(K])-(K]-H)d5di>lin.infi / ([u,]-P^(K])).(P^([u,]-H))d5dt 

+ liminfi / {lu,j~PK{lu,j))-{lu,j-PK{[uej))dSdt>0 (9.4) 

the latter inequality holding due to the properties of the projection operator and the fact that u )^ on Ec . 
Combining (9.3) and (9.4), and rearranging some terms, we readily conclude the weak formulation (4.11) 
of the momentum inclusion. 

Step 2: passage to the Hmit in the semistabiHty condition. It can be performed by the very same 
recovery sequence trick devised in Step 2 of the proof of Theorem 6.1. 

Step 3: passage to the limit in the mechanical and total energy inequalities. It follows from 
(9.2a), (9.2f), (9.2k), and (6.2) that 

$(u(i),z(i)) < liminf $e (Me(t),Ze(i)) for alH e [0,r]. (9.5) 

Combining (9.5) with convergences (9.2) and arguing exactly like in Step 3 of the proof of Theorem 6.1, 
we pass to the limit by lower-semicontinuity in conclude that {u,z,'d) complies for all t G [0,T] with 
the mechanical energy inequality (8.13), with $ in place of $£. Likewise, we conclude the total energy 
inequality (4.12), with $ in place of $£. 

Step 4: mechanical energy equality. Arguing like in Section 8 (cf. [44, Formulae (4.69)-(4.76)]), we 
first of all prove that for all t E [0,T] 

$(u(0, z{t)) - $(uo, zo) + Varx(z; [0, t]) > f (A, u) ds 

Jo 

for any A G L^{0,T;W^-'^{n\Tc;R'^y) with \{t) € du^iu{t), z{t)) for a.a. t e (0,T), (9.6) 

where : W'^^^{n\Tc;R'^) ^ Ty^^^j-^^p^. j^d-j* denotes the subdiffercntial w.r.t. u of the functional 
$ : M^i'2(17\rc;M^) x L°°(rc) M. It follows from definition (4.7) that the operator is given by 

Aea„$(u,z) if and only if Bied^Kiu); W v e W^''^{fl\Tc;M.'^) : 

{X,v) ^ [ Ce{u):eiv) dx + [ Kz[uj-lvj dS + {i,v) , (9.7) 
Jn\rc JTc 
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where we have introduced for notational convenience the functional Ok '■ W^'^{^l\rc',^'^) [0,+oo] 
defined by 3k{u) /x(M) its subdiffcrcntial OOr : VFi'2(^]\rc; R'') ^ W^i-^^^'^p^. j^d^*^ Notice that 
OJk = J* o OIk o J, where J denotes the jump operator J{u) = [u] and J* its adjoint. Now, let us 
observe that 

^ 

'{£,u)ds = 3KHt)) - JkHO)) = lK{lu{t)})^ ^^(["(0)1) =0 

for all e £ L2(0,T;W^i'2(17\rc;M'^)*) such that e{s) G d3K{[u{s)j) for a.a. s e (0,T), 

by the chain rule for the convex functional 3k (cf. [49, Prop. XI. 4. 11]), and by (5.4a) and (4.10). Com- 
bining (9.6)-(9.8), we conclude the following inequality for all t G [0,r] 

$(w(t),z(t)) - $(uo,2o) + Var3^(z; [0,t]) > [ ( [ Ce{u):e{u) dxds + [ [u] • [u] ds) ds . (9.9) 

Jo \ Jn Jtc J 

We now distinguish the two cases ^? = and g > 0. 

Case > 0: arguing by comparison in (4.11) we may readily check that u G L'^{0,T; W^'^{n\Tc;^''^)*), 
which is in duality with u e L^{0,T; W]^'^{il\Tc;M.''-))- As pointed out in Remark 5.3, this entails that ii 
is an admissible test function for the momentum balance inclusion (4.11). Then, upon proceeding with 
such a test and again using (9.8) we conclude for all t G [0, T] that 



f \u{t)\^dx+ f f Be{u):e{u)dxds+ f f Ce{u):e{u) dxds + f f [u] • [u] d^ds 
Jn Jo Jn Jo Jn Jo Jtc 

^ I" 'uQ\^dx + [ ( [ e(i?)B:e(u) dx + / F-udx+ [ f-uds) ds . (9.10) 




2 

" 2 

Combining (9.10) with (9.9), we get the reverse of the mechanical energy inequality, which leads to the 

desired mechanical energy equality. 

Case g ^ 0: From (4.11) one infers that the functional 

e-.v^ [ {Be{u)+Ce{u)-me{d)):e{v)dxdt+ [ kz[u]-Iv] dSdt - [ F-vdxdt- [ f-vdSdt. (9.11) 

Jq Jt.c Jq J^n 

is in L^{0,T;W^^^{n\Tc;R'^)*), and fulfils 

lK{lvj)dSdt> [ lK{luj)dSdt+ [ {e,v^u)dt. (9.12) 

Ec J^c Jo 

Hence, e{t) G d3K{lu{t)j) for almost aU t G (0,T). Thus, (9.8) yields {t,u) ds = for aU t G [0,T], 
which is just relation (9.10) with g ~ Q. Again, we combine the latter with (9.9), and conclude the 
mechanical energy equality. 

Step 5: passage to the limit in the enthalpy equation. It can be developed in the very same way 
as in the proof of Theorem 6.1. We point out that, if (5.8) holds, convergence (9.2i) and (6.8) yield for 
almost all (t, x) G Q the strict positivity of 0(t, x) = Q{i3{t, x)). □ 

Remark 9.1 (Convergence of the reaction force). Notice that, if g = 0, there exists S' > such that, 
for all e > 0, 



= sup 

L2(o,T;VKi.2(n\rc;R<')*) „6L2(o^T;H/i.2(a\rc;a<')) 



(/^)'(M).H dsdt 



which ensues from a comparison in (6.7) by using (9.1) (in spite of the blow-up of the reaction force 
m^yHuS in L°°{0,T;L'^(Tc;R'')), cf. (8.4d)). Passing to the limit in (6.7), we can see that, for the 
subsequence selected in Step 0, J* o (/^)'([[u£]) weakly converges in L'^{0,T;W'^''^{n\Tc;R'^)*) to the 
function i defined in (9.11) . 
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